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THEORETICAL STUDY OF PRODUCTION OF UNIQUE GLASSES IN SPACE 

1.0 INTRODUCTION AND SLMMARY 

This program deals with the glass forming potential of 
unique systems that cannot be obtained by conventional earth 
processing. The objective of this program is to determine if 
the weightless, container less nature of space processing can 
be utilized to produce the glassy form of such materials. 

Glass formation is approached from the kinetic viewpoint. 

The overall objective is to develop analytical functional re- 
lationships describing nucleation and crystal growth in various 
supercooled liquids. The time and temperature dependent rela- 
tionships of nucleation and crystallization (intrinsic properties) 
are used to relate glass forming tendency to extrinsic parameters 
such as cooling rate through computer simulation. 

A thorough understanding of nucleation and crystallization 
kinetics provides a priori knowledge of the ability of a given 
.ystem to form a glass. This will facilitate the development of 
improved glasses by providing a firm theoretical/analytical basis 
for improved manufacturing techniques such as in-space manufacture. 
The ultimate objective of space manufacture is to produce tech- 
nically significant glasses by extending the Earth-limited regions 
of glass formation for certain compositions, or by achieving glass 
formation in other compositions that are not glass formers based 
on empirical Earth observations. 

Single oxide systems were studied initially to aid in devel- 
oping workable kinetic model i and to indicate the primary mater- 
ials parameters affecting glass formation. The literature was 
reviewed and critically analyzed, and kinetic equations were devel- 
oped for homogeneous and heterogeneous nucleation and subsequent 
crystal growth. For homogeneous nucleation the crystal embryos 
are created due to local fluctuations in the structure of the 
liquid phase. If the nucleation is heterogeneous, the phase 
change process is initiated (enhanced) by the presence of foreign 
nuclei. The kinetic relationships derived were applied to known 
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glass formers and non-glass formers. It was found that the 
models qualitatively predict earth-observed behavior for these 
systems . 

i’he theory and analytical expressions developed for simple 
systems then extended to complex oxide systems. The concept 
here ib to develop sufficient analytical relationships to permit 
the prediclicn of the intrinsic glass forming properties of 
unique systems, with the goal of utilizing the in-space environ- 
ment in producr'Ug technically improved glasses. For complex 
systems that da not exhibit glass formation on earth, it is be- 
lieved that the extrinsic effect of heterogeneous nucleation 
often is the major cause. Unfortunately, the analytical treat- 
ment of heterogeneous nucleation phenomena requires much property 
data that are not readily available. To circumvent this problem 
the empirical evidence of the phase change process available 
through North American Rockwell's air suspended laser melting 
experiments (NAS8-28991) was employed. Coupling the information 
present in the filmed melting and subsequent crystallization of 
mullite, the material being studied by Rockwell, with IITRI's 
developed analytics permitted prediction of the intrinsic glass 
forming region of this complex material. The analytical treat- 
ment involved several first order approximations, but the results 
indicate that it is technically feasible to produce mullite 
glasses in the weightless, containerless environment present in 
an earth-orbiting materials processing laboratory. 
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2.0 TECHNICAL DISCUSSION 

The concept of glass forming tendency derives from the 
definition of glass. Morey has defined glass an an "inorganic 
substance in a condition which is continuous with, and analogous 
to, the liquid state of that substance; but which, as a result 
of having been cooled from a fused condition, has attained so 
high a degree of viscosity as to be for all practical purposes 
rigid." Hence, glass-forming materials are ones in which there 
is sufficient transient bonding in the melt to produce a highly 
viscous liquid upon cooling. In the general sense, Morey's 
restriction of glasses to only inorganic materials can be re- 
laxed to include organic materials as well. The supercooled 
amorphous state of aggregation of matter comprising a glass is 
unstable relative to the solid crystalline state. Therefore, 
glass forming tendency is related to the mechanisms and para- 
meters that prevent the liquid-solid transformation from occur- 
ring. Studies of the crystalline transformation can be approached 
from structural, thermodynamic, and kinetic viewpoints. We have 
chosen to adopt a kinetic viewpoint since in general glass forma- 
tion is not related to whether or not a given material can form 
a glass, but rather how fast must the liquid be cooled from its 
fused condition to do so. 

Therefore, we consider a material to be a glass if it can 
be cooled from its liquid state rapidly enough to avoid a certain 
predetermined degree of crystallization. The kinetics of cry- 
stallization of a liquid are determined by two parameters, the 
nucleation rate and the crystal growth rate. The liquid-solid 
transformation occurs by a two-step process of nucleation of 
crystalline embryos and subsequent growth. Nucleation and growth 
rate temperature dependence are illustrated qualitatively in 
Figure 1. Temperature T^^ = T^^^ is the thermodynamic fusion temp- 
erature where the solid and liquid phases are in co-existence. 
Above T^ the material is in the liquid phase. As the liquid is 
supercooled below T^^^, growth can theoretically occur between temp- 
eratures T” and To. However, the embryo nucleation that is 
mo 
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necessary before growth can proceed only occurs between temp- 
eratures T2 and T^. In the temperature region bounded by T2 and 
nucleation and growth occiir simultaneously, i.e., conditions 
are favorable for the complete crystallization transformation pro- 
cess, nucleation and growth. Therefore, the (T2 - T^) tempera- 
ture region is most critical to the prediction of glass formation. 

The present program deals with cha-..acterizing the nucleation 
and growth characteristics of several materials, and presentation 
in a form similar to Figure 1. With the analytical relationships 
describing nucleation and crystal growth, computer -simulated melt- 
queno.h experiments provide information regarding the cooling rates 
necessary for glass formation (i.e., cooling rates necessary to 
pass through the (T2 ~ "^3) temperature region rapidly enough to 
prevent crystallization). 

The kinetic treatment of glass formation necessarily starts 
with derivations of nucleation and growth expressions. In the 
following sections these parameters are discussed in some detail. 
Succeeding sections will deal with the volume fraction trans- 
formed and critical cooling rates, initial computer -mi ^eled re- 
sults on single component systems, a discussion oi the problems 
in predicting the behavior of complex multi-component materials, 
and the use of empirical evidence to permit the analytical treat- 
ment of complex non-glass formers. 
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3.0 KINETICS OF GLASS FORMATION 


The information presented in this section results from a 
critical review of the literature regarding the kinetics of glass 
formation. We have tied together various author's treatments of 
the basic science of the phase change process. This provides the 
basis for our modeling studies on unique systems to follow. 

Equations describing the kinetics of homogeneous nuclea- 
tion, with allowance for transient efl cts, are developed. The 
effect of thp presence of nucleating heterogeneities is also 
discussed. Models for subsequent crystal growth are presented. 
Finally, the method of relating the volume fraction of material 
transformed to the crystalline state to applied cooling rate 
through computer simulation is discussed. 

3 . 1 Theory of Homogeneous Steady State Nucleation 

A supercooled liquid is a raetastable phase relative to the 
solid phase as indicated by free energy considerations. The 
liquid system below the fusion temperature T^^^ tends toward thermo- 
dynamic stability by lowering its free energy through the crystal- 
line transformation. The existence of a supercooled liquid phase 
below Tjjj (i.e., a glass) is the consequence of 1) a surface energy 
harrier between the solid and liquid state, and 2) the kinetically 
inhibited movement of molecules that prevents arrangement in an 
ordered s_, tern (i.e. , crystalline phase) . 

The process initiates as statistical molecular density 
fluctuations causing clustering of molecules or atoms. The 
clusters are called nuclei, embryo particles of solid crystalline 
(transformed) material (transformation to the solid crystalline 
phase can be considered merely as molecular rearrangement into 
an ordered structure) . In tending toward thermodynamic stability 
there is a voltame free energy decrease as an embryo is created 
(i.e., in the transformation from liquid to crystalline). However, 
the formation of an embryo means the formation of a ^oundary, the 
embryo- liquid interface, with a resulting system free energy gain 
due to the interfacial surface energy. Jackson (2) discusses the 
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stability of such a nucleus in an undercooled liquid with 
respect to these two phenomena. For small embryos the surface 
area is large relative to the volume, so that the surface energy 
dominates embryo behavior. Small embryos can decrease the total 
free energy of the system (liquid and nuclei) by shrinking to 
reduce their surface area (i.e., dissolving into the liquid melt). 
Larger size crystalline embryos are controlled by the volume free 
energy term. Large nuclei can reduce the total free energy of 
the system by growing larger, creating more transformed crystal- 
line material (more ordered, more stable). A balance between 
these tendencies defines the critical size nucleus . A nucleated 
embryo smaller than the critical size will dissolve. A nuc’’eated 
embryo larger than the critical size will continue to grow. 

Consider a liquid melt at temperature T with order fluctua- 
tions. This system can be described as a steady state concentration 
of ordered regions (crystalline in structure) of various sizes. 

The change in free energy of the system, AF, due to a local flucua- 
tion creating a spherical nucleus may be expressed: 

AF = - irr^Af^ + Airr^Afg (1) 

where r = embryo radius 

Af^ = free energy difference between the liquid and solid 
state, pci unit voltune 

Af = interfacial (surface) free energy between the phases, 

o 

per unit area. 

The system free energy change has a maximxmi value, AF*, 
for some critical nucleus radius, r*. The critical nucleus 
size, r*, represents the smallest size embryo that can grow with 
a decrease in free energy to form stable nuclei. The critical 
nucleus size is derived as follows. Differentiating AF with 
respect to r gives 

^ = - ^irr^Af^ + STtrAfg (2) 
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Expression (1) has a maximxjm, AF*. for a value of r = r* 
satisfying 


3AF 

3r 


0 


(3) 


Setting the r.h.s. of expression (2) equal to zero and solving 
for r = r* we obtain 


r* 


2Af 

£ 

Tf- 


(4) 


Evaluating the expression for AF (Equation (1)) for r = r* 
yields AF*, the minimum work required to form a stable nucleus: 


AF* 



(5) 


The term AF* is the thermodynamic barrier to the nucleation 
process of forming stable nuclei. We now turn our attention to 
obtaining the surface and volume terms in AF* in mo"^e useful 
form . 


Turnbull (3) has shown that the liquid-crystal surface 
tension, Af^, can be related to the heat of fusion, Ah^, by the 
expression 



Ah^ 

a 


N-1/V2/3 

m 


( 6 ) 


where Ah^ = molar latent heat of fusion of the crystalline phase 
(cal mole ) 

23 -1 

N = Avagadro's number (~ 6 x 10 molecules mole ) 

^m ~ ™olar volxime of the crystalline phase (cc/nlole) . 

The dimensionless terra a relates the proportionality of Ah^ 
and Af , and is constant for a given type of fluid (3). For 

9 

metals a ~ 2, for more complex materials a ~ 3.3. Therefore, 
a has an approximate range 

2 < a < 3.3 (7) 
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The parameter a is defined such that physically the reciprocal 
of a, 1/a, is equal to the nvnnber of monolayers per unit area 
of crystal which would be melted at by an enthalpy AH = Af^. 

The voliome free energy difference between the liquid and 
solid states, Af^, is expressed: 


Af 


r 


= . l^\dT 

j 

m 


( 8 ) 


which is estimated for small supercooling AT = T^ - T: 


Af. 


As^AT 


'm 


(9) 


where As^ is the molecular entropy of fusion (cal mole’^deg"^) 
Since the entropy of fusion is related to the heat of fusion, 
Ahj (cal/mole) 


As , 


Ah^. 
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( 10 ) 
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the expression for Af^ (per unit volume) can be expressed 


Af 


Tm 


(11) 


At large degrees of undercooling AT = T^ - T, this expression 
is modified by the factor T/T 


m 


Af^ (large AT) 


*^m ^m 


( 12 ) 


“•m 


Limiting our analysis to the case of small undercooling, 
use of Equations (6) and (11) for Af and Af„, respectively, 

5 V 

permits expression of r* from Equation (4) in more useful terms 


1/3 


r* = 


2T„V„ 
m m 

aN^/^AT 


(13) 
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For clusters (embryo) below Tjjj with radius r < r*, nuclei form 
and dissolve because the surface energy involved in cluster 
formation is greater than the free energy change accompanying 
solid formation. For r > r*, continuous nuclei growth will occur 
since the surface energy is growing only proportional to r while 
the bulk (volume) free energy terra involved in solid formation 
is growing proportional to r . Above crystalline embryo are 
\jnstable. Critical cluster radii are illustrated in Figure 2 
for various oxide systems. 

The thermodynamic barrier to nucleation, F*, is expressed 

in more useful terms, employing Equation (6) and (11) for Af and 

s 

Af^ respectively: 

,, Ah.T ^ 

AF* = ^ (14) 

^ a-^NAT^ 

It will be convenient to express AF* in terms of another dimen- 
sionless parameter, 6, defined (9) such that 

ASf = NkB = R6 (15) 

where k = Boltzmann's constant 

R = 1.987 cal mole ~^deg~^ = gas constant. 

The heat of fusion, Ah£, thus becomes 

Ah. = As.T„ = R6T„ = Nk3T„ (16) 

r r m m m 

The parameter 3 has range 


1 < 6 < 10 (17) 

with 6=1 for monatomic liquids. More complex structures have 
higher entropies of fusion, with 6 approaching 6 = 10. 

The expression for AF* (Equation lA) is expressed in terms 

of 6: 


AF* 


16 tt 

3 


a^AT^ 


( 18 ) 
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Figure 2 RADII r* OF CRITICAL EMBRYOS AS FUNCTION OF UNDERCOOLING 
-T FOR VARIOUS OXIDES (a = 2) 





Having developed an expression based on thermodynamics and 
free energy considerations for the minimxim work required for homo- 
geneous nucleation (i.e. , the thermodynamic barrier to nucleation) , 
we must now derive an expression for the rate of homogeneous 
nucleation. If it is assumed that critical size clusters (r = r*) 
are formed by statistical molecular fluctuations, the nucleation 
rate will be proportional to a term involving the probability of 
such a fluctuation, exp (-AF*/kT) . The steady state nucleation 
rate involving a Boltzmann- like distribution of critical size 
nuclei is expressed: 

Iq = nv^exp(-AF*/kT) (19) 

where n = number of molecules per unit volume in the system 
N 22 

(n = — ~ 10 ) and v is a molecular attempt frequency for 

''m 

nucleation. 

This result is derived from thermodynamics and thermosta- 
tistics. However, the nucleation process involves molecular 
movement. Therefore, a kinetic or diffusional barrier also 

exists for the nucleation process. The steady state nucleation 

-1 

rate, (nuclei cc sec , is expressed in general terms 

Iq = nv^^exp(-AF*/kT)exp(-AG' /kT) (20) 

where AG' is the activation energy for molecular motion (during 
nucleation) acros.s che embryo-matrix interface. Thus the nuclea- 
tion rate (steady state) is a function of (1) the number of 
molecular or atomistic units available for nucleation, (2) a 
frequency factor describing how often the molecules attempt to 
jump across the liquid-nucleus boundary, (3) a thermodynamic 
barrier to nucleation, the minimum work required to form a stable 
nucleus, and (4) a kinetic barrier to nucleation, involving the 
activation energy for molecular rearrangement. 

If the molecular motion involved in nucleation is treated 
as an activated process (4), a diffusional rate constant for 
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nucleatlon, can be defined such that 

» a^^Vj^exp(-AG'/kT) (21) 

where a^ is the molecular jximp distance. 

Combining Eqtiations (20) and (21) yields the following 
expression for the steady state nucleation rate: 

nD 

Iq = — I exp(-AF*/kT) (22) 

o * 


The problem now arises as how to evaluate the nucleation 
diffusional rate constant, D^. For liquids which crystallize 
without change in composition (single component systems, for 
instance) long range diffusion processes are not required. All 
of the atomic or molecular units required for the ordered crystal- 
line structure are in the local vicinity of the liquid-crystal 
interface. This type of transformation is termed non-reconstruc- 
tive. The activation energy, AG' , for diffusion at the liquid- 
crystal interface will be roughly the same order of magnitude as 
the activation energy for viscous flow. This is the case for a 
liquid transforming to solid without change in composition since 
the movements of the atoms or structural units on the nucleation 
surface is similar to the reorientation of structural \inits and 
bond switching in the flow of a viscous liquid. Therefore, in 
this simple case, the nucleation diffusional rate constant, D^, 
V7ill be approximately equal to the self diffusion coefficient of 
the undercooled liquid: 


n s 


(23) 


The liquid self diffusion coefficient, D , is related to the 
bulk viscosity, r\, through the Stokes -Einstein relation 


kT 


STia^n 


( 24 ) 


where is the diameter of the diffusing species 


T 
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Therefore, for the case of a single component substance, 
or a more complex liquid that crystallizes with no change in 
composition, the steady state nucleation rate is expressed (for 
small degrees of undercooling below T^) 


I 


o 


nkT 

STra^^n 


exp 


- 16 itST„^l 

m 

3a%Al2 . 


(25) 


where all terms have been previously defined. 

The general temperature dependence of the homogeneous 
nucleation rate described by Equation (25) is shown in Figure 1. 
For small degrees of undercooling below T^, AF*, the thermo- 
dynamic barrier, is large since Af^, the volume free energy 
change in the transformation is small. This results in low 
nucleation rates. With further supercooling Af^ increases 
until AF* ^ AG' , the kinetic barrier to nucleation. This condi- 
tion results in maximum nucleation rate. For a large AT the 
nucleation rate decreases to a negligible level as AG'> > AF*. 

In the case where a compositional change accompanies the 
liquid- solid transformation, the steady state nucleation rate 
(Equation (25)) must be modified to account for the long range 
diffusion processes that are required for molecular rearrangement. 
This subject will be treated in Section 5.0. 

3.1.1 Time -Dependent Homogeneous Nucleation 

In this discussion transient effects on nucleation are 
-.onsidered. Glass forming tendency is related to how fast the 
system can pass through the region of simultaneous nucleation 
and growth (i.e., the (T 2 - T^) region shown in Figure 1). The 
time spent at any particular temperature level may then be less 
than the time required to establish a steady nucleation rate 
(i.e., the time required to build up the required Boltzmann -like 
distribution of embryos) . 
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Incorporation of time -dependence into nucleation theory 
has been treated by Hillig (5) and Hammel (4) , and takes the 
form 


I 


t 



(26) 


where * transient nucleation rate 

I * steady state nucleation rate 
t = time . 

H'itnmel (4) has expressed the parameter t (after Collins (6)) 
as : 


144kTT V 
m m 

N^/^aa ^as.AT^D 
o r 


(27) 


where D is the appropriate diffusion coefficient (given by the 
Stokes-Einstein equation for transformation without change of 
composition) and all other terms have been previously defined. 
This equation was derived for the ideal case of instantaneous 
cooling from to T. 

Hillig (5) derived an approximate expression for x using 
random walk diffusion theory to describe the mean time t to build 
a critical nucleus: 

u 2 *2 

TTV, r* 

T = t = - ^ X - . - (28) 

4DV 

m 


where and are the mole volumes of the liquid and precipi- 
tating phase, respectively, D is again the appropriate diffusion 
coefficient (D = D for crystallization without change in composi- 
tion) , and X is the mole fraction of the precipitating phase in 
the liquid (i.e. , X = 1 for crystallization without change in 
composition) . 

Hillig' s t parameter has been shown (4) to more accurately 
describe time-dependence of the nucleation rate, based on 
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correlation with experimental results. Therefore, we shall 
employ Equations (26) and (28) when simulating nucleation in 
computer generated melt-quench experiments. 

The nature of transient effects on nucleation rate is 
illustrated qualiuatively in Figure 3. Due to diffusional effects 
successively longer times are require-^ to achieve steady state 
nucleation at successively lower temperatures. If the hypotheti- 
cal material whose nucleation characteristics are shown in 
Figure 3 is cooled from T^ to T^ detectable nucleation would be 
avoided at any cooling rate*. However, a detectable level of 
homogeneous nucleation occurs between temperatures Tj^ and T 2 . In 
order to avoid this nucleation in a melt-quenching experiment the 
temperature region T^^ - T^ must be passed through in less than 
time t 2 - If this glass were successfully quenched to temperature 
Tg, for instance, it could be held at this temperature level for 
a period of time tg before detectable nucleation would occur 
(e.g. , for annealing). 

As shall be discussed in Section 5.0, the transient 
nucleation formulation provides one possible way of dealing with 
multicomponent systems where long range diffusion processes are 
required for molecular rearrangement. 

3.2 Heterogeneous Nucleation 

Thus far we have only treated the case of homogeneous 
nucleation in pure substances free of impurities , insoluble 
particles, etc. Foreign surfaces present in the liquid, such 
as container walls or insoluble particles, however, tend to 
reduce the barrier to nucleation represented by the surface 
energy between the liquid and solid phases. As a result the 
critical embryo size is reduced and the "supercooling temperature" 
is raised. The "supercooling temperature" is the temperature T 2 
where the first detectable nucleation is ob'^erved upon cooling 
from the melt. This effect is shown qualitatively in Figure 4. 
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At heterogeneous surface sites the computation of the 
critical embryo size becomes complicated because of the various 
surface energies involved and the particular geometry of the 
embryo surface. Furthermore, the wettability of the foreign 
surface by the liquid must be taken into account. 

For our present purposes we shall only treat the effect 
of heterogeneous nucleation qualitatively, following the treat- 
ment of several authors (references 2, 4, 5, 7-9) in order to 
illustrate the nature of the analytical expressions. 

The work term for heterogeneous nucleation, , 

given generally as 


AF*^ - f(0) • AF* 


(29) 


where the factor f(6) is expressed 


f(9) 


(2 cos9)(l - cos6)^ 
4 


(30) 


The angle 6 is a function of the balance of surface tensions and 
is given by: 


COS0 



(31) 


where a j , o and o, , are the interfacial energies between the 

S Li s p L p 

solid-liquid, solid-impurity surface, and the liquid-impurity 
surface, respectively. 

The lowering of the net interfacial energy in heterogeneous 
nucleation depends on how well the impurity surface wets the 
nucleating phase in the presence of the liquid. 

Therefore, if heterogeneous nucleation sites are present 
the minimum work required ( thermo djmamic barrier to nucleation) 
to form a staole nucleus is reduced. The treatment is further 
complicated by the fact that the concentration of nuclei deriving 
from a heterogeneous mechanism will be some function of the concen- 
tration of impurity sites present, and not given simply by a 
nucleation rate formulation. 
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3. 3 Crystal Growth Kinetics 

Once formed by homogeneous ox heterogeneous nucleation 
mechanism, stable nuclei will continue to grow at a rate determined 
primarily by the rate at which the atoms necessary for growth can 
diffuse to the surface of the crystal, and by the ease with which 
they can free themselves from the attractions of their neighbors 
in the liquid phase and form new bonds in the specific positions 
determined by the structure of the growing crystal (10) . The 
fundamental concept here is that growth ts considered in terms of 
molecular rearrangement. 

Crystal growth is either interface-controlled or diffusion- 
controlled. Interface-controlled growth occurs when the rate 
controlling step occurs at the liquid-crystal interface. This 
normally occurs for single component glasses, for instance, 'here 
only short range molecular rearrangement at the liquid-crystal 
interface is necessary. Diffusion-controlled growth occurs when 
the rate controlling step is the long range diffusion of a given 
species in the bulk liquid. Thus, diffusion-controlled growth 
would normally apply in the case where crystallization is accom- 
panied by p. large change in composition, such as for complex 
multicomponent systems. 

3.3.1 Interface-Controlled Growth Kinetics 

For this work we shall restrict ourselves to the case of 
crystallization without change in composition as we have done 
for the case of homogeneous nucleation in the preceding sections . 
Thus, we will treat in this section only the mechanisms and 
kinetics of interface-controlled growth. 

Two basic mechanisms have been proposed for interface- 
controlled growth; continuous growth and lateral growth (11). 

In continuous growth the crystal interface advances by molecular 
incorporation which can occur ,ith equal probability everywhere 
(except for certain anisotropic effects) . Lateral growth occurs 
either by a two-dimensional nucleation mecharism or by a screw 
dislocation spiral growth ramp mechanism depending on the crystal 
perfection. 
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In systems which crystallize without change in composition 
the nature of the liquid-crystal interface strongly influences 
the kinetics and morphology of crystallization. Different models 
for interface controlled growth are each based on a different 
assvunption concerning the interface and the nature of the sites 
on the interface where atoms are added or removed. Again, it is 
important to think of crystal growth as merely a molecular 
rearrangement process. 

The general form of the growth velocity for all interface- 
controlled growth processes is: 

y = fa^Vgexp(-AG'7RT)[l - exp(-AG/RT) (32) 

where f = fraction of preferred growth sites on the interface 
(i.e., fraction of sites available for growth, 

0 < f 1 1). 

V = frequency factor for molecular transport at the 
liquid-crystal interface (during growth) . 

a^ = distance advanced by the interface in a unit kinetic 
process (~ molecular diameter) . 

AG = free energy change accompanying the liquid-crystal 
transformation 

AG” = free energy of activation or kinetic barrier for 
the movement of an atom across the liquid-crystal 
interface during crystallization (i.e., growth). 

The kinetic term AG” involved in growth is not necessarily 
equal to (or even the same order of magnitude) as the kinetic term 
involved in nucleation, AG' (Equation (20)). Growth may be 
governed by atomic movements from a great distance away from the 
liquid-crystal interface. Nucleation (initiation of the trans- 
formation) is governed by atomic movement relatively close to the 
developing nucleus. 

Assuming that we can trea*- the molecular movements involved 
in crystal growth as simply activated processes, we can define a 
diffusion coefficient, or rate constant, for the molecular 
rearrangement 








D = a ^v exp(-AG'VRT) (33) 

o ^ o 

in a manner similar to the case of nucleation described in 
Section 3.1. The general interface-controlled growth expression 
(Equation (32)) thus becomes: 


[l - exp(~AG/RT)] 


Experience on a variety of pure substances which crystal- 
lize without change in composition has indicated that the slow, 
rate controlling step in the interface growth process is the same 
molecular process that occurs in the bulk liquid above the thermo- 
dynamic fusion temperature, Tjj^(3) . Therefore, in this most 
straightforward of cases, the rate constant for interface- 
controlled growth, D , will have approximately the same value as 
the liquid self diffusion coefficient, D , and thus be approximated 

O 

by the Stokes-Einstein relation: 


= D 


kT 

3traQn 


The free energy change accompanying the liquid- crystal 
transformation, AG, is given by 


AG = ASjAT 


The growth velocity for interface- controlled growth without 
composition change is thus expressed 




r ,-ASfAT| 

1 - exp — — — I 


or in terms of the B parameter (As^ = NkB = 


since 


y = 


STTa^^n i 


I'-BATlI 

“ii 


R = Nk 
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We will now briefly treat the various models that have 
been proposed for interface-controlled growth. The models and 
discussion presented have been taken from Uhlmann (12) . 

3. 3. 1.1 Normal (Continuous) Growth 

In normal or continuous growth atoms can attach to or be 
removed from any site on the interface. Thus, there are no 
preferred growth sites in the interface, and f in Equation (37) 
becomes unity. 


y 


kT 

3Tra^^n 


1 


exp 


; -As^AT; 

i ± — 

\ RT I 


(40) 


For small departures from equilibrium, this model predicts a 
linear relation between growth rate, y, and undercooling, AT. 

For this model to correlate with experimental data the liquid- 
crystal interface must be rough on an atomic scale. 

3.3. 1.2 Screw Dislocation Growth 

In the screw dislocation model growth occurs at step sites 
provided by screw dislocations intersecting the interface. The 
fraction of preferred growth sites, f in Equation (37), is 
expressed 


f 


aoAhfAT 


4TraT V 
m m 


2ttT, 


m 


(41) 


where it has been assumed that only molecular transport within 
a molecular diameter, a^, of the dislocation ledge results in 
attachment. For small undercooling, AT, this model predicts a 
growth rate which varies with AT^. For the screw dislocation 
model to apply, the interface must be smooth on an atomic scale, 
and be relatively imperfect in the crystallographic sense. 

3. 3. 1.3 Surface Nucleation Growth 

According to this model growth takes place at step sites 
provided by two-dimensional nuclei formed on the interface. The 
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growth rate can be expressed 


y = Avexp (-B/TAT) 


where the exponential corstant is 

T 

B = o IP m E 


kAh- 


( 42 : 


(43) 


The term o^. is the specific edge surface energy of the nucleus . 
The frequency factor, v, can be derived from the growth rate 
constant : 



The growth velocity predicted by this model should vary exponen- 
tially with the undercooling, AT, and for small AT should be 
unobservably low. For this model to correspond to experimental 
observation the interface must be smooth on an atomic scale and 
be defect-free. 

This comment and the foregoing discussions regarding 
expected morphology for each of interface-controlled growth 
models will have application in the experimental stages of our 
NASA program (to be discussed later) . For materials which have 
not been extensively studied, morphological studies (SEM, etc.) 
will indicate, by inference, which growth model applies to a 
given material. Utilizing empirical results of melt-quench 
experiments to correlate with our kinetic relationships will 
help to indicate how space processing can be utilized to produce 
glasses unobtainable on Earth. 

3.3.2 Significance of Entropy of Fusion 

As discussed by Uhlmann (12) growth rates depend critically 
on the molecular structure of the liquid-crystal interface. The 
interface structure, in turn, depends significantly on a bulk 
thermodynamic property, the molecular entropy of fusion, As^. 
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For materials characterized by low entropies of fusion 
(ASj < 2R) the liquid-crystal interface should be rough on an 
atomic scale, defects should not affect growth, and normal growth 
kinetics are predicted (i.e. , Equation (40)). On the scale of 
light microscopy the crystal- liquid interface should be non- 
faceted. 

For materials characterized by large entropies of fusion 
(As£ > 4R) faceted interface morphologies should be observed, 
defects should be important to growth, and the kinetics of the 
normal growth model should not apply. Except as limiting cases, 
the obs-?.r-^ed growth velocities for large As^ materials should 
not agree well with behavior predicted by screw dislocation 
model or the surface nucleation model either. 

3 . 4 Volvime Fraction Transformed and Critical Cooling Rates 

Glass forming tendency can quantitatively be described 
in terms of the volxrae fraction of crystalline material formed 
during a certain quenching operation. Uhlmann (13) has expressed 
the volume fraction, v^, crystallized in time, t, (for small v^) 
as 


Vj - l/3TTly^t^ (45) 

where I = nucleation rate (nuclei cc ^sec"^) 
y = growth velocity (cm sec”^) 
t = time (sec) . 

This relation is valid for interface-controlled growth, and thus 
applies to single component or congruently melting (i.e., without 
change in composition) compounds. 

The critical cooling rate for glass formation, i.e., the 
minimum cooling rate necessary to avoid a certain volume fraction 
of transformed material, can be estimated by a procedure employed 
by Uhlmann (13) . T-T-T (time-temperature-transformation) curves 
are constructed as illustrated qualitatively in Figure 5. Using 
Equation (45) , the time required to transform a given amount of 
material is calculated as a function of temperature. A volxmie 
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Figure 5 TIME-TEMPERATURE-TRANSFORMATION CURVE 
(TAKEN FROM UHLMANN (13)) 
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fraction transformed * 10 ^ is regarded as just-detectable. 
The data are presented graphically as AT vs t in the T-T-T 
plot. The critical cooling rage can be approximated (13) by 
the expression 

dT = (46) 

3t . Y“ 

critical 

where AT„ = T„ - T„ 

N m N 

Tj^ = temperature at the nose of the T-T-T curve 
T», = time at the nose of the T-T-T curve. 
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4.0 COMPUTER MODELED RESULTS FOR SELECTED SINGLE 
CO MTotffiNT SYSTEMS 

The kinetics developed for nucleation and growth are now 
applied to several materials systems. The well-characterized 
silica, Si 02 , system is investigated for the purpose of deter- 
mining the relative importance of the various kinetic and 
thermodynamic parameters in controlling transformation behavior. 
The ^ 2^3 system is studied to determine if the derived kinetic 
relationships will predict its good glass forming qualities , and 
imply a physical reason for the observed behavior. The alumina, 
AI 2 O 2 , system is studied to determine if the kinetics can predict 
the poor glass forming tendency that is observed on earth, and 
imply a physical reason for this phenomenon. 

4.1 The Silica System (Si02) 

The pertinent materials properties for the silica system 
to be used in the derived k'.netic nucleation and growth expres- 
sions are tabulated in Table I. Viscosity data for tho silica 
system are shown in Figure 6, compared with viscosity data for 
other oxide systems. 

The materials parameters were inserted into the derived 
nucleation and growth kinetics (Equations (25) , (26) , (28) , 
and (37)), with results illustrated in Figures 1 and 8. The low 
temperature cut-off in the nucleation ra<-e has bean discussed 
previously. Figure 7 indicates the relative time scale of this 
effect. A nucleation rate of 1 nuclei cc~^sec~^ has arbitrarily 
(but by popular custom (5)) been taken as the level of detectable 
nucleation. As illustrated in Figure 7, detectable homogeneous 
nucleation for SIO 2 does not occur upon cooling from the melt 
(T^) until roughly 1700° to 1750°K. Above about 1200°K the 
thermodynamic barrier to nucleation controls behavior. Below 

~ 1200°K the kinetic (diffusional) factor dominates . The peak 

0 

growth rate for Si02 (Figure 8) is roughly 40 A/min, and occurs 
at roughly 50°K undercooling belcw T^. The crystal growth rate 
is seen to be zero at T^, the thermodynamic fusion temperature. 
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TABLE I 


MATERIALS PARAMETERS FOR SiO, 


Ahj 

m 

B 

f 

X 


2000 cal mole 


2000°K 


-1 


Ah 

As, 


^ = leal mole'^K"^ 


.5 


2.5 A ( ~ 0-0 distance in SiO^ tetrahedron) 
1 (normal growth assximed) 

1 (pure substance) 


= V- = 27.6 cc mole 


-1 


n (viscosity) : see Figure 6 

a =2.5 (representative value) 

23 -1 

N = 6 X 10 molecules mole 

N 97 -1 

n = — = 2 X 10 molecules cc 

^m 

AT = T„ - T 
m 

k = Boltzmann's constant = 1.38 x. 10”^^ erg C 

3.3 X 10'^^ cal C’^ 

R = gas constant = 1.987 cal mole"^deg~^ 



Temperature (°C) 
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For small values of supercooling (temperatures during quench 
below Tjjj) , the free energy term is dominant and increases with 
incieasing supercooling (i.e., the growth rate increases with 
decreasing temperature) . At large degrees of supercooling 
below the kinetic factor (diffusion) begins to con‘ ol, and 
the growth rate begins to decrease with decreasing temperature. 

At temperatures far below the diffusion kinetics dominate: 
the extremely high viscosity inhibits crystal growth, and the 
growth rate drops to zero. 

The crit5cal temperature region for glass forming tendency, 
the region of simultaneous nucleation and growth, is illustrated 
qualitatively in Figure 9 in a "T 2 - T^" region plot. If homo- 
geneous nucleation only is considered, it is observed that Si 02 
glass formation depends on the ability to pass through the temper- 
ature range ~ 1600 - 1700°K rapidly. The transient behavioi 

shown in Figure 7 indicates that this critical region must be 

3 

tranversed in a time less than approximately 5 x 10 sec to avoid 
the liquid-crystal phase transformation. This time scale limit 
is readily attained by commercial earth processing methods. 

4.2 Hypothetical Glass-Parametric Study 

The computer software that has been developed to generate 
nucleation and growth behavior using the kinetic and thermodynamic 
relationships derived in Section 3.0, provides a convenient means 
of determining which material parameters are of primary importance 

in determining glass forming tendency. This information is 
critical to our overall task of investigating the glass forming 
jndency of unique systems, or in synthesizing systems for space 
manufacture . 

For this purpose a hypothetical glass was po.stulated that 
possessed the viscosity- temperature behavior of Si 02 . and the 
parameters a and 6 were varied over realistic ranges (i.e., 

2 ± cx ^ 3, I 561 10). The resulting parametric graphical 
representations are presented in Figures 10 and 11, and serve to 
indicate the importance of these parameters on the (steady state) 
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Figure 9 T,-T. REGION NUCLEaTION AND GROWTH PLOT FOR SILICA 




Figure 10 STEADY STATE NU CLEAT ION RATE FOR A HYPOTHETICAL GLASS (T^ - 2000"K, 
P - 1, a - VARIABLE); THERMODYNAMIC AND KINETIC BARRIERS 




Figure 11 STEADY STATE NUCLEATION RATE FOR A HYPOTHETICAL GLASS (Tm - aOOO'K, 
a - 2, B - VARUBLE); THERMODYNAMIC AND KINETIC BARRIERS 
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nucleation rate. The horizontal line 3 !:he level of detectability 
of nucleation (assumed to be 1 nucleu per cm per sec) . The 
temperature at which the nucleation curves intersect the hori- 
zontal level of detectability line is temperature T2 shown in 
Figure 1. Above T2 no nucleation should be observable. Since 
nucleation and growth occur simultaneously between T2 and T^, 
we are interested in decreasing (T2 - T^) ; therefore, decreasing 
T2 will enhance glass forming tendencies. Additionally, we are 
interested in material parameters that tend to decrease the slope 
of the nucleation rate -temperature relationship in the region 
just below temperature T2 (this is important due to limitations 
in cooling rates attainable) . 



Figure 10 illustrates the effect of a for the hypothetical 
glass with the B parameter fixed at unity. It is observed that 
decreasing a decreases T2, and thus tends to increase glass forming 
tendency . 

Figure 11 illustrates the effect of B for the hypothetical 
glass, with the a parameter fixed at a value of a * 2. It is 
observed that the B parameter has a larger effect on the nuclea- 
tion rate than does the a parameter. Temperature T2 decreases 
with increasing B • 

It is apparent that for enhancing glass forming tendencies 
large B and small a parameters are required. This is based on 
the observed decrease in the temperature T2 and decrease in the 
slope of the nucleation curve near T2. This can also be seen by 
considering the expression for AF*, the free energy of a critical 
embryo, which represents the magnitude of the barrier to nucleation. 
A large B represents a material with a large heat of fusion; a 
complicated structure with a large amount of energy involved in 
the phase change. A small a represents a material with a large 
heat of fusion relative to the solid- liquid interfacial free 
energy. Values for B (or As£) for several systems are shown in 
Table II. 
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TABLE II 

PERTINENT PROPERTIES Of VARIOUS MATERIALS 


Material 

As^ 

(cal/mole °K) 

Ah^ 

(cal/mole) 

6 

T„(°K) 

Si02 

1 

2D00 

.5 

2000 

®2°3 

8.2 

5900 

4.1 

723 

GeO 2 





a-Al«0« 

11.3 

26000 

5.6 

2300 


Lithium Silicates 5 

Potassium Silicates 12 

CaO-Al 203*2 Si02 16.1 

Soda-Lime Glass 8 

Metals 2.3 
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Thus far we have only discussed thermodynamic parameters 
a and B. Since there also exists a kinetic barrier to nuclea- 
tlon and growth, it would be expected that viscosity, n, would be 

a parameter of primary importance. Uhlmann (13) has discussed 
the importance of two kinetic parameters on the tendency of a 
material to be a glass former; (1) a high viscosity at the 
fusion temperature, and (2) a rapidly rising viscosity with 
decreasing temperature below (i.e., large dn/dT). The latter 
is related to the position of T^^^ relative to the glass transition 
temperature, T . We will not elaborate on the parameters n and 
dn/dT in a parametric study. However, they will be considered 
in subsequent discussions of the glass forming tendencies of 
B2O2 and AI2O2. 

4.3 The B2O2 System 

In this section nucleation and growth kinetics are applied 
to the B2O2 system to predict the observed good glass forming 
quality for this material, and to imply a physical reason for this 
phenomena. Qualitative indications of glass forming tendency 
can be obtained by considering only steady state nucleation, 
ignoring for the moment transient effects. Steady state data 
represent an upper bound for nucleation behavior. 

The pertinent material properties for the 820^ system are 
tabulated in Table III. The viscosity of ^ 2 ^^^ illustrated in 
Figure 12. Over the temperature range 250-500°C the literature 
data shown were averaged and fit to a 6th degree polynomial 
expression: 


log n = AT^ + BT^ + CT^ + DT^ + ET^ + FT + G 

(47) 


with resulting coefficients: 
A = 2.24 X 10"^^ 

B = -4.56 X 10“^^ 

C = 3.82 X 10"® 

D = -1.73 X 10"^ 
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TABLE III 

MATERIALS PARAMETERS FOR B 2 O 3 



Ah^ 

= 

5900 

cal mole”^ (reference 21 ) 

h 

T^rr, 


450°C 

= 723°K 

TCI 




T 

j 

ASc 


Ahc 

= 8.2 cal mole'^K"^ 

1 . 

f 


*^111 


1 " 



As 


1 . 

B 

— 

r 

R 

= 4.1 

}■ 


= 

0 

2.5 A 

(assvimed) 

M ft 

f 

= 

1 (assumed normal growth) 

1 ' 
i. 

X 

= 

1 (pure substance) 


a: variable, 2 < a < 3 


r 

r 

« > 
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Bruckner (22) 

Macedo (23) 
Napolitano, et.al. (24) 
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E =■ 4.81 X 10"^ 

F = -8.9 X 10"^ 

G = 9.77 X 10^ 

The computer-generated viscosity data for 820^ are also shown 
in Figure 12 for comparison. 

These material properties were inserted into our computer 
software and the steady state nucleation rate for 820^ computed 
and plotted for variable a as illustrated in Figure 13. The 
computed growth rate for 82©^ is illustrated in Figure 14. Glass 
forming tendency for 82©^ is deduced from considerations of the 
nucleation and growth behavior shown in these two figures. 

Referring to the growth curve (Figure 14) , it is observed 
that the growth rate for 82O3 is zero at T^, peaks at roughly 

700°K, and decays to zero by roughly 625°K. The steady state 
nucleation rate for 82©^ (Figure 13) was computed for the known 
B parameter (B = 4.1) and for a probable range of tne a parameter. 
It is observed that for all values of a, the nucleation rate is 
always below the detectable limit (assumed 1 nucleus per cc per 
sec), and thus effectively zero. Therefore, although 82©^ 
exhibits a theoretical growth rate much higher than Si02 (see 
Figure 8) , critical size crystal embryo are never homogeneously 
nucleated. These analytical results predict, therefore, that 
82O2 would be an excellent glass former since it can not be made 
to nucleate and crystallize when cooling from the melt, regardless 
of the cooling rate. This prediction is supported by empirical 
evidence indicating that crystallization of 820^ from a dry melt 
has never been observed (12), 

The parameters of primary importance in describing glass 
forming tendency from a kinetic standpoint are 1) the magnitude 
of the viscosity at T^^^, 2) the slope of the viscosity- temperature 
function below T^, dri/dT, and 3) the magnitude of the molecular 
entropy of fusion As^ (related to AH^ and B) . The former two are 
related to the diffusional barrier and the latter to the thermo- 
dynamic barrier to the phase change process, liquid state to 
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Figure 13 STEADY STATE NUCLEATION RATE FOR B,0, FOR VARIOUS VALUES OF THE 



Figure 14 CRYSTAL GROWTH RATE F(» B,0. 
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crystalline state. The magnitude of the melting point viscosity 
for B2O2 is orders of magnitude less than that of Si02 (see 
Figure 6) . This would indicate that nucleation-crystallization 
should occur more readily for B2O2 (lower diffusional, kinetic 
barrier). Since this is not the case, the factors inhibiting 
the phase change in B2O2 are the higher molecular entropy of 
fusion (8.2 for B2O2, roughly 1.0 for Si02> , and the steepness 
of the viscosity- temperature relation below T^^^ for B2O2 (dn/dT 
greater for B2O2 than for Si02 as shown in Figure 6) . 

To illustrate the relative importance of As^ and dn/dT 
in inhibiting nucleation and crystallization in the B2O2 system, 
we have used our computer system to describe the nucleation 
behavior of a hypothetical liquid with the viscosity- temperature 
relation of B2O2, but with the molecular entropy or SiO,, (i.e., 

As£ = 1 ) . Our hypothetical glass then exhibits the kinetics of 
a material that will not nucleate and crystallize (B2O3) and the 
thermodynamics of a material that will nucleate and crystallize 
(Si02) . The nucleation behavior of this hypothetical material 
is shown in Figure 15, with the behavior of B2O2 (As£ = 8.2) 
also shown for reference. It is observed that the nucleation 
rate of the hypothetical liquid is now above the detectable limit. 
Thus by lowering the thermodynamic barrier to the phase change 
for B2O2 (i.e., As£ or B) we are able to produce a hypothetical 
glass that will nucleate. The growth rate for this hypothetical 
liquid was computed and is shown in Figure 16. Comparing the 
relative positions of the nucleation and growth curves for this 
hypothetical liquid (Figures 15 and 16) indicates its glass 
forming qualities. For convenience these nucleation and growth 
curves are qualitatively sketched in a "T2 - T^ region" plot in 
Figure 17. It is observed that a small region of simultaneous 
nucleation and growth exists between 625 and 650°K. Thus, crystal- 
lized B2O2 could be obtained upon cooling from the melt (by 
holding the liquid between 625 and 650°K) if the molecular entropy 
of fusion (or heat of fusion) could be reduced substantially. 

This hypothetical example is presented as an indication of how we 
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Figure I') STEADY STATE Nl CLEAflON RATE FOR B2O3 AND A HV JTHETICAL LIQUID (as^ gj^Q 






Figure 17 NUCLEATION AND GROWTH CHARACTERISTICS OF A HYPOTHETIC iL 
LIQUID (Dg Q , ^siOo’ 



might eventually be able to synthesize materials to exhibit 
desired nucleation-crystallization tendencies. 

4.4 The AI2O2 System 

Nucleation and growth kinetics also have been applied to 
the AI2O2 system. This system exhibits very poor (non-existent) 
glass forming tendencies by conventional earth methods. The aim 
here is to predict this poor glass forming tendency of AI2O2 
ar.d imply a physical reason for this phenomenon. 

The pertinent materials properties for the alumina system 
are tabulated in Table IV. The major difficulty in describing 
nucleation and growth kinetics for the Al20^ system is the complete 
lack of viscosity data below the fusion temperature, T^. This 
situation exists since AI2O2 never been observed on earth 
as a glass (i.e., at temperatures below T^^^) . Above the melting 
point, however, data have been reported ( 17 ), and are illustrated 
in Figure 6. With no reference point such as the glass transi- 
tion, T , available, the problem remains as to the shape of the 

. O 

viscosity- temperature curve for Al20^ below T^^^. For our purposes 
we have chosen two cases; 1) extrapolating below T^^^ with the 
shape of the 62^3 (i.e., steep dn/dT) and 2 ) extrapolating 

below Tjjj with the shape of the Si02 curve. These two cases thus 
represent the upper and lower bounds for dn/dT of the various 
oxide: shown in Figure 6. For both cases the extrapolated curves 
were expressed algebraically using curve fitting techniques for 
use in our computer software. 

The homogeneous nucleation and growth characteristics 

of AI2O2 are summarized in the "T2 - T^ region" pi t shown in 

Figure 18 . Behavior depicted by solid lines represents the case 

wheie AloO- is assumed to have a dn/dT behavior below T similar 
z j m 

to that of Si02 (see Figure 6) . The dotted lines in Figure 18 
represent the case where the viscosity-temperature dependence 
of B2O2 is assumed below the AI2O2 fusion temperature, T^^^. Since 
AI2O2 has such a low viscosity at its fusion temperature (melting 
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TABLE IV 

MATERIALS PARAMETERS FOR AI 2 O 3 


l! 


? * • 


I . 
I ■ 


m 

f 

Ah, 


2300°K 


= 2.5 A (as Slimed) 


= 1 (assumed) 


-1 


= 26,000 cal mole (reference 21 ) 


1 * 
i I 


-1 -1 

LSf = = 11.3 cal mole K 

^ T„ 
m 


As 


= 5.65 


1 ! 

I I 

t 

\ 




a 


= 2 , 2.5, 3 (variable) 
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f 
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Figure 18 HOMOGENEOUS NUCLEATION AND GROWTH CHARACTERISTICS OF AI 2 O 3 


1 


point) , theoretical growth rates are much higher than for either 
pure SiO« or Both of the Al^O- growth rate curves shown 

Z ^ J ^ J 1 Q (J 

in Figure 18 have peak magnitudes of ~10 A/nin, compared to 
40 A/min for Si02 (Figure 8) and 2 x 10^ A/min for B20^ (Figure 14) . 
To illustrate the very large theoretical growth rate of AloO- , the 

AO ^ J 

peak growth of B^O^ (2 x 10 A/min) is reached by the alumina 
system within - 10"^ undercooling below T . The 40 A/min peak 

Ul Q 

growth rate of Si 02 is reached by AI 2 O 2 within ~ 10~° under- 
cooling below T . 

m 

At large degrees of undercooling below T^^^, in the diffusion- 
controlled portion of the growth curve, the steepness of the 
viscosity-temperature relation is seen to determine the tempera- 
ture at which the growth rate drops to zero. Since we are dealing 
with tbe approximate upper and lower boundaries of the possible 
viscosity relation of alumina below T^^^, this temperature (T^ in a 
conventional "T 2 - T 2 region" nucleation and growth plot) lies 
between approximately 1300° and 2050°K. 

The steady state homogeneous nucleation rate for AI 2 O 2 is 
also shown in Figure 18 (for a variable) , For the condition where 
the AI 2 O 2 viscosity has the general temperature dependence of 
B 2 O 2 below T^ (steep dn/dT) it is observed that homogeneous nuclea- 
tion is below the detectable limit. Where AI 2 O 2 is assiamed to 
behave more like Si 02 , which is more probable, it is observed that 
a detectable level of homogeneous nucleation occurs for a = 3. 

The shaded region shown in Figure 18 thus represents the region 
of simultaneous (homogeneous) nucleation and crystal growth for 
AI 2 O 2 system. 

If homogeneous nucleation and subsequent growth are con- 
sidered, we have shown that to avoid crystallization upon cooling 
from the melt the temperature region 1300° to 1625°K must be 
passed through rapidly. In practice, however, AI 2 O 2 is known to 
crystallize almost immediately upon cooling below the fusion 
temperature, T^ (2300°K) . The reason for this discrepancy is that 
our analytical predictions of glass forming tendency have not 
accounted for heterogeneous nucleation. In the case of hetero- 
geneous nucleation, any insoluble impurity or external surface 
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will serve to lower the size of a critical embryo and thus 
effectively increase temperature T 2 (the temperature where 
detectable nucleation first appears upon cooling from the melt) . 
This effect is shown in Figure 19. Temperature T 2 , corresponding 
to homogeneous nucleation, is shifted to a higher temperature 
(T^') if heterogeneous nucleation is possible (i.e., if extrinsic 
nucleation sites are present). As illustrated in Figure 19, 
simultaneous heterogeneous nucleation and growth can (qualita- 
tively) occur at temperatures just below the fusion temperature, 

Tjj^. Most investigators of crystallization phenomena indicate 
thru nucleation appears to initiate heterogeneously in most 
materials. Thus an extrinsic property controls the glass forming 
tendency of AI 2 O 2 . at least on earth where the complete elimina- 
tion of heterogeneous nucleation sites is not possible. Complete 
elimination of external nucleation sites in Al 20 ^ is necessary 
if AI 2 O 2 glass is to be obtained since the growth rates are so 
high. Therefore, AI 2 O 2 might be a good candidate for space 
manufacture since processing could be performed container less , 
with no external surfaces acting as nucleation sites. Before 
the space-processing candidacy of AI 2 O 2 is determined, however, 
several areas must be investigated in more detail. These include: 
1) glass quality regarding crystal size and concentration if 
homogeneous nucleation only is possible (i.e., considering the 
high growth rates and attainable quench rates from 1625*^ to 
1300°K, perhaps enough crystalline phase would be nucleated 
homogeneously to yield a poor quality glass even with space 
processing) , and 2) insoluble impurity levels attainable in 
AI 2 O 2 precursor materials (i.e., perhaps enough impurity sites 
will be available for heterogeneous nucleation to make elimina- 
tion of surface nucleation sites (crucible wall) through space 
processing only a second order improvement) . 

Clearly, if space processing is to be employed to eliminate 
heterogeneous nucleation sites leading to high quality glasses, 
then the mechanisms of heterogeneous nucleation must be well 
known before candidate materials can be chosen with any confidence. 
This presents a very formidable problem since a purely analytical 
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treatment of heterogeneous nucleation requires much information, 
such as the various surface energies discussed in Section 3.2, 
that is not readily available. Section 6.0 discusses an empirical 
approach that has been used to circumvent this problem. 

4.5 Summary of Modeling Results for Single Component Systems 

The kinetic and thermodynamic relationships describing 
glass formation that were developed in Section 3.0 have been 
applied to three relatively simple systems: 1) Si 02 , a well- 

characterized system, 2) B 2 '’ an excellent glass former, and 
3) A12021 a material exhibiting no glass forming tendency based 
on empirical earth observation. It was found that viscosity, 
the viscosity- temperature relation, and the entropy of fusion 
are the parameters of prime importance in determining glass 
forming ability. For Si02 and £ 20 ^, our analytical treatment 
predicted glass forming behavior that is consistent with the 
known behavior of these materials. For the AI 2 O 2 system it was 
determined that the extrinsic effect of heterogeneous nucleation 
must be investigated in order to accurately predict the level of 
improvement in glass forming ability that could be obtained 
through in- space processing. 

These results gave us confidence that we could now proceed 
to more complex systems covering more unique materials, and 
indicated to us the critical areas of investigation for such 
materials . 
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5.0 APPLICATION OF DERIVED KINETICS TO MORE COMPLEX SYSTEMS 

The nucleation and growth kinetics that have been developed 
in Section 3.0 apply in general to single component substances 
or congruently melting (i.e., without composition change) com- 
pounds. The liquid-solid transformation process for these 
materials is termed non-reconstructive . No interatomic bonds 
within the participating molecules need be broken, and only short 
range diffusion processes are required for the transformation to 
occur at the liquid-crysti*! interface. In this instance, the 
molecular movements are treated as simply activated processes, 
with rate constants approximately equal to the coefficient of 
self-diffusion in the bulk liquid. This is convenient since the 
self-diffusion coefficient is related to the viscosity through 
the Stokes-Einstein equation, and viscosity data are more readily 
attainable kinetic data than diffusion data. 


In multicomponent systems, however, the liquid to solid 
transformation often involves bond breaking and/or long-range 
diffusion processes, as discussed previously. In this instance 
the tra.'S format ion is termed reconstructive. In network liquids, 
such as Si 02 » interatomic bonds must be broken in the network 
prior to molecular rearrangement. Since this bond breaking in 
network liquids must also precede viscous flow or self-diffusion, 
the free energy of activation for these processes will also be 
applicable as the kinetic barrier to the phase change process (25) . 
Thus, in th. reconstructive transformation of a network liquid we 
can approximate the activation energy and the rate constant in the 
same manner as for non-reconstructive crystallization: 

AG' (nucleation) = AG" (growth) = AG (flow 

dL 

activation energy) ,, 

{HO) 

and, 


D 


n 



D 

s 


kT 

STia^n 

o 


(49) 
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However, in reconstructive crystallization where a large 
change in composition is involved, long range diffusion processes 
are required to bring the appropriate atomic species to the 
liquid-crystal interface. In this case the rate limiting step may 
be the diffusion of a particular (i.e. , the slowest moving) species 
in the liquid matrix. In dealing with nucleation and growth 
kinetics and subsequent glass forming tendency of such materials , 
several methods have been postulated. Uhlmann and Chalmers (8) 
suggest that for nucleation involving large changes in composi- 
tion that the kinetic barrier to nucleation, AG' in Equation (20), 
should be taken as the activation energy for diffusion of the 
slowest moving component in the matrix, and that the pre-exponen- 
tial factor, nv, should be reduced by the mole fraction of the 
precipitating component. This presents a problem since the 
appropriate diffusion data are not as readily available as bulk 
liquid viscosity data. Furthermore, the kinetic term for nuclea- 
tion, AG' , is not necessarily equal to or even the same order of 
magnitude as the kinetic term for growth, AG" in Equation (32). 

This is due to the fact that growth is governed by atomic movements 

from a great distance from the interface, whereas the molecular 
movement involved in the nucleation process occurs relatively 
near the liquid-crystal interface. In general, growth in multi- 
component systems where composition is a variable tends to be 
diffusion-controlled rather than interface-controlled. Hammel (4) 
has discussed the formulation for computing the volume fraction 
of material transformed for diffusion controlled growth: 


_ 8tt q3^.3/2^ .5/2 


(50) 


where D' is the diffusion rate constant and S is a supersatura- 
tion term. Again, however, diffusion coefficients are required 
that are not readily available. 


The question arises, then, of how we can predict the 
glass forming tendency in systems where diffusion controlled 
behavior is expected due to large compositional changes during 
the transformation, and for which the required diffusion data 
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are not available. Hillig (11, 26) has proposed that incorporation 
of a transient nucleation rate (51 will adequately account for 
long range diffusional effects that occur in a reconstructive 
transformation. Hammel (4) has applied transient nucleation 
analysis (Equations (26), (27), (23)) and diffusion controlled 
growth analysis (Equation 50)) to predict the volume fraction of 
cristobalite (Si02) precipitating from E glass (13Na02 ' llCaO- 
76Si02) during cooling from the melt. In the nucleation rate 
expression Hammel used a value of X = .1 (mole fraction of pre- 
cipitating phase in the melt. Equation (28)) to account for the 
long range molecular rearrangement. A value of X equal to unity 
is employed for a pure single component substance or a congruently 
melting compound. The determination of the appropriate value for 
X to use in a given reconstructive transformation is not clearly 
defined, however. Hammel (27) and Hillig (26) have proposed that 
X be determined by considerations of 1) what the expected race- 
limiting species will be and 2) its concentration in the melt. 

In Hammel 's treatment of diffusion controlled growth described 
above, a value of S in Equation (50) was determined using 
Frank's (28) tabulated values for diffusion controlled growth. 

Unfortunately, the state-of-the-art in this nucleation- 
crystallization area has not reached the point where predicted 
transformation kinetics correlate well with reality in all 
instances. Especially in the case where large composition 
changes accompany the transformation, it may be that empirically 
derived nucleation and growth kinetics can be employed as dis- 
cussed by Uhlmann (12) . We have had the opportunity to investigate 
this. Section 6.0 describes how empirical evidence can be employed 
to describe glass forming tendency in a case where sufficient 
data is lacking to permit a purely analytical approach. 
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6.0 ANALYSIS OF COMPLEX SYSTEMS USING EMPIRICAL EVIDENCfi 

The homogeneous nucleation and growth kinetics that have 
been treated thus far apply generally to single component sub- 
stances or congruently melting compounds. The last section 
discussed several possible methods of treating less ideal mater- 
ials which \mdergo a composition change upon crystallization. 
However, even for simple substances we have found that a major 
problem exists in the manner in which we can handle the effects 
of heterogeneous nucleation. Referring to Section 3.2, we have 
shown that the analytical treatment of heterogeneous nucleation 
requires information about the magnitude of the solid-liquid, 
solid-impurity, and the liquid-impurity interfacial energies. 

This information is not readily available. Thus, analytical 
determination of the intrinsic glass forming ability and space 
candidacy of a material such as AI 2 O 2 which is believed to cry- 
stallize via a heterogeneous nucleation mechanism, is a formid- 
able task (See Section 4.4, The AI 2 O 2 System). 

For complex systems, i.e., systems for which our ana-.ytics 
are not entirely adequate, we now must rely, in part, on empirical 
evidence of the phase transformation process. Tne follov/ing 
section details how we have employed the empirical evidence avail- 
able in a film of the crystallization of mullite to separate out 
heterogeneous nucleation, and determine the region of intrinsic 
glass formation. The empirical evidence mainly consisted of an 
indication of the viscosity- temperature relationship during the 
phase change process. Mullite was chosen for this analysis since 
to date it is the only material that such a film is available for. 
Ideal mullite (3Al202*2Si02) crystallizes without composition 
change, so we can apply the same general kinetics that we have for 
single component systems. 

It is our eventual aim to develop this analytical-empirical 
approach to the stage where it can be employed on even more complex 
materials. Once the mechanisms controlling the Earth-limited 
regions of glass formation in such unique materials have been de- 
termined, we can more readily assess how a technique such as 
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in-space processing can extend the boundaries of glass forma- 
tion yielding technically improved products. 

6 . 1 The Mullite System QAl pO^ ‘ ZSiO p) 

Mullite glass has only been obtained (on earth) in very 
small quantities by splat cooling techniques. Our eventual aim 
is to determine a) the intrinsic glass forming region of mullite, 
b) the causes of the earth-observed poor glass forming ability 
of mullite, and c) whether or not a space-processing technique 
could be used to produce large mullite optical elements . 

The glass forming tendency of the mullite ( 3 AI 2 O 2 • 2Si02> 
system is indicated by the relative positions of the nucleation 
and growth kinetrc', as discussed previously (i.e., T 2 -T 2 region 
of simultaneous nucleation and growth). Unfortunately, an analy- 
tical treatment of nucleation and growth cannot be initiated 
until we have some basic kinetic and thermodynamic data for this 
system; mainly the heat of fusion (ov the entropy of fusion) and 
the viscosity below the fusion temperature, Tm. An extensive 
literature search turned up neither the required entropy of 
fusion data, nor the required viscosity data. The only kinetic 
information we have for this system is 1) the fusion temperature, 
and 2) the viscosity at the fusion temperature. We do, however, 
have some empirical evidence of the crystallization of mullite: 
a filmed record of the air suspended laser melting and subsequent 
crystallization of mullite. This empirical evidence of the 
phase change process gave us information regarding the temperature- 
viscosity relation where the process initiated, and led to the 
determination of the entropy of fusion of mullite. The information 
in the mul?tite crystallization film therefore plays an integral 
part in our determination of the glass forming behavior of mullite, 
as will become evident as vre proceed through our analysis. 


* Performed by R. A. Happe, North American Rockwell Corp . , under 
NASA Contract NAS8-28991. 
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6.1.1 Determination of Critical Property Data 


Our analysis starts with an estimation of the viscosity 
of glassy mullite. We have employed the same procedure we used 
when investigating the AI2O2 system; extrapolation 01 _he vis- 
cosity below the liquidus temperature (i.e., Tm) by assuming the 

shape of the known SiO_ n-T curve. Viscosity at the liquidus 

^ (29'i 

temperature was taken from Takamori and Roy'* ^ . The resultant 
viscosity relation is shown in Figure 20. It is emphasized here 
that this represents a first order approximation of the mullite 
viscosity. No attempt was made to consider arguments regarding 
the structural similarity between mullite and silica. Such more 
rigorous treatments of the viscosity of mullite below T based 
on available structure-viscosity theories will be considered later. 


To proceed with our analysis we must now determine the 
entropy of fusion of mullite. This parameter was shown to be of 
prime importance in determining glass forming behavior. To derive 
this information "2 first need to have an estimate of the crystal 
growth rate of mullite. R. A. Happe analyzed fi?ms of the 
laser melting of mullite and presented overall growth rate data^^^^ 
for each experimental melt. The overall crystal growth rate (ad- 
vancement of crystallization front) averaged for all experimental 
runs was 


y '\i 


776 — -V 6x10^ 
sec 


A 

min 


For the purposes of our analysis we shall assume that this repre- 
sents a peak growth rate, Up. (Refer to the shape of the y-T 
relationship shown in Figure 8, for instance). This assumption 
is believed to be a valid first order approximation based on the 

following argument. Our previous work has shown that the peak 

10 ° 

growth rate of AI2O2 is « 4x10 A/min. The liquidus temperature 

viscosity of alumina ('\^ 0.66 poise, reference 31) is a factor of 

3.18 lower than that of mullite (2.1 poise, reference 29). It 

v^ould thus be expected that the crystal growth rate of AI2O2 would 

be greater than that of mullite by roughly this factor. Multiplying 
90 10 o 

the 6x10 A/min mullite growth rate by 3.18 gives 'v- 2x10 A/min 

for a rough estimation of the AI2O2 growth rate. This estimated 
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growth rate for AI 2 O 2 , based on the mullite films, is within 
a factor of two of the AI 2 O 2 growth rate we have computed from 
analytical considerations. Thus, the mean growth rate observed 
in the mullite laser melting films is quite close to the value 
for the peak growth rate. 

Continuing our analysis, we now have a peak growth rate 
for mullite, and must now arrive at the temperature below 
where mullite is expected to exhibit its maximum growth rate. 

A rough indication of this temperature is obtained by compar- 
ison with other systems. This comparison is based on observ- 
ations made of the growth rate temperature dependence of the 
Si02, ^2^3’ previously studied. In order 

to approximate growth rate temperature dependence, three 

temperatures are defined: T , T , and T . T is the thermo- 

m p o m 

djmamic fusion temperature, is the temperature where the 

growth rate peaks at a meiximum value, and T^ is the low-side 

temperature where the growth rate decays to a negligible 

magnitude. It was found that definite correlations exist 

for the temperature ratios T /T , T /T , and T /T for each 

o m o p pm 

of the systems studied, as shown in Table V. The data for 
these computations vere extracted from the computed growth 
rate temperature curves for Si 02 , ^2^3’ ^^2^3 

previously. At uming that this correlation is valid for the 
mullite system, we now compute the temperature where mullite 
is expected to exhibit maximum growth: 

Tp = .97(T^) = .97(2123°K) = 2059°K (51) 

The viscosity that glassy mullite would exhibit at this 
temperature (see Figure 20) is obtained graphically as 
n - 6.3 Poise at 2059°K. 

It is emphasized here that we are not aware of any 
theoretical basis for this correlation method for obtaining 
the approximate temperature where mullite will evhibit maximxim 
crystal growth. In view of the '' ack of sufficient analytical 
or experimental data, however, .. ,s empirical/analogical 
approach is the only one that will allow us to go forward 
with the analysis. Our overall approach to describing the 
glass forming behavior of complex systems is an iterative one. 
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TABLE V 



CORRELATIONS OF CRYSTAL GROWTH BEHAVIOR 
FOR THE Ai 20 ^, B 2 O 3 . and SiC 2 SYSTEMS 


>roximate Temperatures and Ratios (°K) 


Material 

T 

0 

"p 

m 

p m 

0 m 


Si 02 

1600 

1940 

2000 

.97 

.8 

.83 

®2°3 

620 

700 

725 

.97 

.85 

00 

AloO. 

Z j 

1980 

2230 

2300 

.97 

.86 

.89 


I , ; ~ thermodynamic fusion temperature 

I 

f Tp = temperature of peak growth rate 

I ■ ~ low-side temperature where the growth 

• , rate is 'v zero (or y < < y peak, 

M i.e., the tail of the curve) 





If our final results are not compatible with experience and 
judgment, we go back etrd examine the credibility of the 
assumptions that have been made along the way. 

Proceeding with our analysis, we can now evalute our 
analytical growth rate expression at the peak growth temper- 
ature (2059°K) to compute the only other unknown parameter 
for mullite, the molecular entropy of fusion. Assiuning 
normal interface-controlled growth kinetics (3 to 2 mullite 
melts without composition change) . the crystal growth rate 
expression is (see Section 3. 3. 1.1); 


fkT 



1 


exp 


^-As^AT ' 
RT~ 


(52) 


The pertinent parameters are: 

u = .776 (growth rate from mullite film) 

T = 2123°K 
m 

T = 2059°K (Temp, where p = y peak) 

f = 1 (congruent melting) 

o 

a^ = 2.5A (assumed) 

Ti = 6.3 poise at 2059°K (see Figure 20) = 

15.06x10® 

cm. 

R = 1.987 cal/mole°K (gas constant) 

AT = T -T = 64°K 
m 

k = 3.3x10 cal/°K (Boltzman's constant) 

These values were substituted into Equation (52) to evaluate 
the entropy of fusion, As^, for mullite: 

As. = 6.78 (53) 

^ mole°K 

The heat of fusion of mullite is calculated to be; 

Ahf = =. 14,400 11^ (54) 


65 



These computed thermodynamic parameters for mullite agree 
well with values derived from estimates of the heat of fu ion and 
entropy of fusion based on weighted averages of the mole constitu- 
ents (references 32 and 33) : 

^(Al 20 ^) + ^(Si 02 ) = mullite 


•^^^^f/Al 203 ) -^^^^f(Si02) 

.6(26,000) + .4(2000) = 16,400 
or Ah^ (mullite) = 16,400 


and thus 


A s^- (mullite) 


16,400 cal 

Tte" ■ moTe^ 


Thus, the value of the entropy of fusion of mullite we have com- 
puted from the empirical results of laser melting experiments is 
within roughly 107o of the value computed from a first order ap- 
proximation. 

6.1,2 Nucleation and Growth Kinetics 


We have thus far determined an approximate viscosity-temp- 
erature relationship for vitreous mullite, and utilized the films 
of the laser melting/crystallization of mullite to derive an ap- 
proximate value for the entropy of fusion. With this tnermodynamic 
and kinetic information we can now proceed to analytically des- 
cribe the glass forming tendency of mullite. This is ■'e by 
obtaining a rough idea of the region of simultaneous nu .'eation 
(homogeneous) and crystal growth for this system. 


For our initial analysis of this system we need only investi- 
gate steady state nucleation behavior. We have shown previously 
that for the congruent melting of the 3Al202-2Si02 system, the 
steady state homogeneous nucleation rate is given by the expres- 


sion: 


I = 
o 


nkT 

3-rra^^ri 


exp 


16 

-16it f m 

— 3- — \ 7 

^ Ra' TAT 


(55) 
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For the mullite system the following materials parameters apply: 


As£ = 6.77 cal mole"^ K”^ 

a = 2 , 2 . 5 , 3 (variable) 

T = 2125°K 

m o 

a^ = 2.5 A (assximed) 
n(T)= See Figure 20 
n,k,R = Constants 


These values were inserted into Equation (55) and the steady 
state homogeneous nucleation rate for mullite computed and 
plotted as illustrated in Figure 21 (for variable a) . 


The crystal growth rate for mullite (normal interface con- 
trolled growth being assumed) is given by the expression 


y = 


fkT 

STra^ n 


1 - exp 


■As^AT 

~^T 


(56) 


The pertinent material parameters (given above, with f equal to 
unity) were inserted in this relation, with the resultant growth 
velocity-temperature curve shown in Figure 22. The extremely 
large theoretical growth rates shown for mullite (up to 6x10^ 
X/min) are the result of the low liquidus temperature viscosity 
(see Figure 20) . 


Superimposing the homogeneous nucleation and growth rate 
behavior of mullite on che same temperature scale will indicate 
the "T 2 -T 2 " region of simultaneous nucleation and growth (i.e. , 
a qualitative indication of intrinsic glass forming tendency) . 

This is illustrated in Figure 23, where for convenience the growth 
velocity is expressed as log y. It i s observed that if the a 
parameter for mullite lies somewhere in the range 2.5 £ a £ 3 
(reasonable values) , then it would be expected that mullite would 
intrinsically nucleate and crystallize upon cooling from the melt 
at a temperature between 'v 1500°K and 1700°K, and would continue 
to undergo the liquid-solid transformation down to a temperature 
a- 1200°K. Thus, in the absence of nucleating heterogeneities 
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Figure 21 . Nucleation Rate for Mullite as a Function of Alpha 
(Steady State) 






Figure 23. Intrinsic Nucleation and Growth Characteristics for Mullite 

(Shaded region ■" simultaneous nucleation and growth) 



the temperature region of roughly 1700°K-1200°K would have to 
be passed through rapidly to form glassy mullite (i.e,, avoid 
detectable crystallization) . 

In practice, however, it is a common laboratory occurrence 
that mullite crystallizes immediately upon passing the liquidus 
temperature (thermodynamic fusion temperature, T^) . The crystal- 
lization of mullite can be routinely observed at undercooling 
temperatures (i.e., degrees below Tj^^) much less than the 400°K 
undercooling predicted in Figure 23 (i.e. , Tj^-T 2 = 400°K) for 
homogeneous nucleation (i.e., intrinsic glass formation). This 
is the result of nucleating heterogeneities which tend to shift 
the nucleation curve, i.e., T 2 , closer to T^ = T^j. Heterogeneous 
nucleation is the process by which the liquid-solid phase trans- 
formation is initiated by foreign surfaces present in the liquid 
such as container wall asperities or insoluble particles. These 
surfaces tend to reduce the barrier to nucleation represented by 
the surface energy between the liquid and solid phases. As a 
result of these extrinsic parameters the critical embryo size is 
reduced and the supercooling temperature is raised (the super- 
cooling temperature i.s the temperature where the first detectable 
nucleation is observe . upon cooling from the melt) . 

This discrepancy between our analytical results and general 
experience for the mullite system was also observed in our analy- 
sis of AI 2 O 2 (refer to Section 4.4). However, for mullite the 
laser melting films provide a unique way of dealing with hetero- 
geneous nucleation in the absence of the required analytical data. 
This will be discussed in the following section. 

6.1.3 Treatment of Heterogeneous Nucleation 

The air suspended/ laser melting experiments provide a con- 
venient means of investigating the heterogeneous nucleation effect. 
In the actual laser melting of mullite performed by R. A. Happe, 
several phenomena were observed or recorded. We have reviewed 
films of this experiment at IITRI and have observed two important 
phenomena. First, crystallization of the mullite sphere pro- 
ceeded from a finite number of nucleation sites (from one to three) . 
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Second, the nucleation and initial crystallization appeared 
to occur on the surface of the mullite sphere. This was dis- 
cernible as a contoured crystallized surface observed while 
the mullite sphere was rotating. These phenomena indicate that 
the nucleation was heterogeneous, not homogeneous. 

The following argument supports our contention that since 
only one to three sites were observed, the nucleation had to 
have been caused by a heterogeneous mechanism. From the ob- 
served rate of advancement of the crystallization front, .V cm/ 
sec) and the sample size ('^ .6 cm diameter), the total time 
elapsed for the process was on the order of 0.5 sec. The cooling 
rate for the experiment was estimated by Happe to be 100°C/sec. 
Thus, during the observed transformation the temperature decreased 
roughly 50°C. Assuming a uniform temperature distribution and 
referring to the mullite nucleation rate shown in Figure 21. it 

is illustrated that if the nucleation process were homogeneous, 

2 3 

then roughly 10 -10 sites would have nucleated. During the 
process only 1 to 3 sites were observed, supporting the hetero- 
geneous nucleation hypothesis. 

R. A. Happe 's analysis of the mullite films and the actual 
melting experiments provide additional information^ ’ ' about 

the crystallization behavior of mullite; 1) an undetermined 
amount of initial superheating above T^^was employed in all ex- 
periments, 2) a temperature drop of 600*^C to 1000°C was observed 
from the initial superheated condition to the point where crystal- 
lizatiei was first visually detected, and 3) before and during 
crystallization the viscosity of the sample appeared low and did 
not rapidly change with decreasing temperature. The latter item 
was inferred from the observed rapid oscillating motion, "jiggling,” 
of the mullite sphere. 

Using this information about the freezing of mullite, we will, 
attempt to describe the experiments shown in the mullite films 
analytically. The idea here is to see how intrinsic the glass 
formation was with the aim of eventually being able to predict 
how we can conduct an experiment in space that will yield a glass 
that cannot be produced on Earth. 
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One of the most basic parameters of glass forming ability 
is the temperature T 2 where the first detectable phase change 
is observed upon cooling from the melt. If the crystallization 
of mullite observed in the films were produced by intrinsic 
homogeneous nucleation, then the temperature at which the 
first detectable phase change would be detected (upon cooling 
from the melt) would be 1700°K (T 2 in Figure 23). Figure 20 
indicates that the viscosity of amorphous mullite at 1700°K 
would be 1x10^ Poise. This level of viscosity is believed 
to be too high to produce the "jiggling” observed in the 
mullite films. A viscosity of 10^ Poise is inconsistent , we 
believe, with R. A. Happe's observations that at the time of 
initial crystallization, the viscosity was very low, and did 
not increase significantly with decreasing temperature. A 
viscosity of lO^Poise is three orders of magnitude greater 
than the liquidus temperature viscosity. Therefore, the 
viscosity of the mullite in the films just prior to crystall- 
ization must have been much lower than 10^ Poise, say, in the 
range of ^ 10^ Poise. This means that the temperature of 
initial detectable crystallization would have to have been 
roughly 1900°K (from Figure 20, n = 10^ Poise at 1700^^0 . 

Thus temperature T 2 in the actual melting experiment would 
be 1900°K, rather than 1700°K as predicted by Figure 23 for 
the case of homogeneous nucleation. 

Shifting the initial crystallization temperature T 2 
upwards to 1900°K shows the effect of het rogeneous nucleation. 

This is illustrated in Figure 24. The nucleation curve is 
shifted from T 2 to T 2 ' , closer to the fusion temperature 
T^ = Tj^. This evidence for heterogeneous nucleation cor- 
relates with our observations that nucleation appeared to in- 
itiate at a finite number of locations on the mullite surface. 

If, then, the initial crystallization occurred at approximately 
190C°K (T 2 '), the amount of undercooling becomes T^^^ - T 2 ' “ 200°K. 
Therefore, the amount of superheating in these experiments 
would be in the range of 400'^ to 800°K above the liquidus 
temperature T^ (v2100°K) . 
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This is based on R. A. Happe's observations that the mullite 
temperature dro,ned from 600° to 1000°C from the time the 
laser was shut off until crystallization was first detected 
visually. 

6.1.4 Critical Cooling Rates 

For the intrinsic case of homogeneous nucleation, our 
first-order approximation analysis has r.hown that when slow 
cooled from the melt (T^ = =2125°K) one would expect to 

observe the liquid-solid phase transormation when the 
temperature reached the 1500° to 1700°K level, depending on 
whether the a parameter is closer to 2.5 or 3.0 (see range 
of T 2 indicated in Figure 23) . Tht laser melting films were 
interpreted to estimate tnat for the actual crystallization 
observed by Happe, the first detectable phase, transformation 
occurred at roughly 1900°K. This illustrates the effect 01 
heterogeneous nucleation sites, anc shifts the nucleation 
curve upwards from T 2 to T 2 ’ as shown in Figure 24. 

We now attempt to quantify these cases by computing the 

cooling rate required to avoid detectable phase transformation. 

(13) 

Using Uhlmann’s procedure ^ , which was discussed in Section 

3.4, computer software V7as developed to perform this analysis. 

A typical computer-generated time-temperature-transformation 
curve is illustrated in Figure 25. 

The first set of computations of critical cooling rates 
assumed that a = 2.5, as shown in Figure 26. The critical 
cooling rate for glass formation was computed for three cases: 

(A) for homogeneous nucleation, i.e., intrinsic behavior, 

(B) for nucleation that is heterogeneous to the extent that 

we have previously estimated the laser melting experiments to 
be (i.e., T 2 shifted upwards to T 2 ' ® 1900°K) , and (C) for 
nucleation that is heterogeneous to the extent that might be 
present in conventional earth melting experiments (i.e , T 2 
shifted upwards to T 2 ' 2100°K, which corresponds to nucleation 

just below the fusion temperature (2125°K) from a site such 

as a crucible wall ox an undis solved refractory impurity 
particle). These three cases are illustrated in the ’'T 2 -T 2 " 
region plot shown in Figure 26. Figure 27 illustrates these 
three cases repeated for the case where a * 3.0. 
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Extrixislc and Intrinsic 

Figure 26. Nucleation and Growth Characteristics for Mullite <a«2.5) 
(Shaded region ■ simultaneous nucleation and growth) 






(shaded region ■ simultaneo”*) nucleatlon and growth) 




The results of this analysis are presented in Table VI. 

The 10^-10^°K/sec critical cooling rate predicted for a "con- 
ventional" Eai-th-bound experiment seems reasonable since "splat" 
techniques have been found to be the only methods capable of produc- 
ing significant glassy mullite phases on Earth. Quench rates 

C £ 

of 10'^-10”OK/sec are estimated for splat cooling methods. 

The cooling rate of 10'^-10^°K/sec that we have predicted 
to be critical for the conditions under which the laser melting 
of mullite was performed also appears reasonable. This is evi- 
denced by the fact thai glassy mullite was not produced in these 
experiments, where the experimental cooling rate was estimated to 
be only v 10^°K/sec^^^^ . 

The computed critical cooling rate for mullite free of 

1 3 

nucleating heterogeneities is 10 -10 °K/sec, as shown in Table 
VI. Conventional quenching methods such as radiation cooling lie 
within rhi*" region. Therefore, according to this analysis, mul- 
lite has intrinsic glass-forming qualities realizable through 
relatively conventional quenching methods. Achieving intrinsic 
glass formation thus becomes a problem of eliminating hetero- 
geneous nucleation sites. This can be partially accomplished by 
going to containerless space processing where no sample-crucible 
interface exists that would provide heterogeneous sites. 

6.1.5 Summary of Mullite Analysis 

To circumvent the problem of the lack ol kinetic data for 
the complex mullite system, the empirical evidence of the phase 
trc^ns format ion process available in the North America! Rockwell 
films was utilized. Our conclusions regarding the laser melting 
experiments and the glass forming ability of mulliit , derived 
from our "To-T^" region nucleation and growth analysis are as 
follows ; 

• the liquid- solid transformation observed in the films 
initiated by a heterogeneous nucleation mechanism. This 
was deduce.! "rom the fact that crystallization initiated 
on the surface of the sample sphere from a finite number 
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of sites. If the nucleation were homogeneous, the entire 
sphere volume would have nucleated simultaneously. 

• the first detectable crystallization occurred at a 
temperature level, T 2 ' of 1900°K (i.e., at an under- 
cooling of 200°K) . This was deduced from the observa- 
tions of R. A. Happe that at the time of initial vis- 
ually detectable crystallization, the mullite sphere ex- 
hibited a "jiggling" motion. Such motion would be pro- 
duced by a material with a viscosity of 10^ poise, 
corresponding to a temperature of 1900°K for mullite. 

• an initial superheating of 'v- 400°K existed in the laser 
melting experiments. This was computed from Happe 's 
estimate of the degree of cooling obtained after laser 
shut down to the first detectable phase change, together 
with the estimated temperature T 2 '. 

• the cooling rate for the laser melting of mullite was 

2 

V 10 °K/sec (derived from the cooling curves shown by 
Happe) . 

• for the amount of heterogeneous nucleation existing in 
the laser melting experiment, a cooling rate of 10 to 
10^°K/sec would be necessary for glass formation. This 
result was computed from IITRI's analysis, and is con- 
sistent with the estimated cooling rates of these exper- 
iments. and the fact that mullite glass was not obtained. 

• for the amount of heterogeneous nucleation present if a 
crucible were used or if no superheating were used, a 
cooling rate of v 10^°K/sec would be required for glass 
formation. This result was computed from IITRI's analy- 
sis, and is consistent with the known cooling rates of 
splat cooling techniques (v lO^-lO^oK/sec) , nnd the fact 

that significant glassy mullite phases have only been 
obtained using splat techniques, 
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• a cooling rate of 10 to 10 OK/sec will be necessary 
for glass formation in the mullite system if all hetero- 
geneous nucleation sites are eliminated (represents the 
intrinsic behavior of mullite, computed from IITRI's analy- 
sis) . The cooling rate required for intrinsic glass forma- 
tion in mullite is within the range of relatively conven- 
tional quenching techniques such as radiation cooling. 

This first order approximation analysis has indicated that 
if a technique such as weightless, containerless space processing 
could be utilized to remove all heterogeneous nucleation sites 
from mullite, then mullite glasses which are relatively unobtainable 
on earth could be produced. The potential exists for relatively 
large optical elements since conventional cooling methods could be 
employed. 

It is emphasized here that this analysis has given us a first 
order approximation .f cooling rates critical for glass formation. 
Enough evidence has been ■'resented, we believe, to warrant further 
consideration of mullite as a space processing candidate material. 
However, several refinements to our analytical technique should 
eventually be considered. First, transient nucleation has not 
been considered in this analysis. We have only included I^, the 
steady state nucleation rate. Secondly, when we shifted the 
nucleation curve upwards to simulate the existence of heterogeneous 
nucleation sites, the shape of the nucleation- temperature curve 
was not changed. This mainly involves not accounting for changes 
in the shape of the viscosity-temperature curve at different temp- 
erature levels. Finally, considerations should be given to theor- 
etical viscosity-structure relationships to improve upon the extra- 
polation technique used to estimate mullite viscosity below the 
fusion temperature. Nonetheless, the present analysis has exhibited 
sufficient correlation with empirical observations to be of value 
in predicting the range of cooling rates necessary for intrinsic 
glass formation. Further discussion of possible refinements to our 
analysis is presented in Section 7.0. 
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6 . 2 Additronal Experiments - Calclmn Alumlnate Glasses 


In addition to the mulllte work described above, various 
calcium aluminate compositions were investigated. High alumina 
calcium aluminate has possible ir-transmission applications, but 
it is not a good glass former based on earth-boxand experience. 

It is our eventual aim to determine if space processing could be 
utilized to extend the regions of glass formation for this material. 

The xCa 0 *yAl 20 ^ compositions for this study were prepared 
by reaction sintering of a cold pressed precursor. The precursor 
materials were calcium carbonate and aluminum tri-hydrate. The 
general reaction equation is : 

xCaC 03 + y [a 1203 ' 3 H 20 ] xCa 0 -yAl 203 + xC02f+ 3yH20^ (57) 

The phase diagram for this system and the earth -determined glass 
forming region are illustrated in Figure 28. Particular composi- 
tions 'or our study lie within the glass forming region and on the 
high-aiumina side of the glass forming region - as designated by 
# 1 , # 2 , etc. 

Peferring to Figure 28, compositions #1 and #2 (i. 3 ., 50w/o 
alumi-'.a and 70 w/o alumina, respectively) were prepared at IITRI 
and then air suspended and laser-melted at Rockwell Corp, The good 
glass forming tendency of a 50 w/o AloO- composition was confirmed 
in the laser melting experiment^ . Additionally, the 70 w/o 
alumina composition was also found to exhibit glass formation^^^\ 
This is significant since this composition lies outside the pre- 
viously accepted glass forming region shown in Figure 28. 

Tlie scheduling of the preparation of the 75 w/o and 80 w/o 
alumina samples at IITRI was not compatible with the time schedules 
of the Rockwell program. In the future, we intend to prepare the 
75 w/o and 80 w/o alumina systems. After sintering at roughly 
1600”c the existence of the desired AI 2 O 3 + 3 Ca 0 - 5 Al 203 material 
(See Figure 29) will be verified by X-ray diffraction. The sample 
will then be sent to R. A. Happe for subsequent laser melting. 




Liqoidus temperature (*C) 


Compositions Studied 


1 2 34 

I 111 



. The Glass Forming Region in Relation 
to the Phase Diagram for the Calcium 
Aluminate System (38) 
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Eventually, the alumina content of the IITRI -supplied calcium 
aluminate samples will become high enough that good quality glass 
will not be obtained in the laser melting experiments. At that 
point, we can apply our analytical nucleation and crystallization 
relationships to the films of the crystallization phenomena as 
we have done for the mullite system. When this work is completed 
the results will be reported to NASA in the form of supplementary 
information. In this manner, the intrinsic glass forming ability 
can be determined, and the candidacy of the particul.:.r calcium 
aluminate composition can be assessed. 



7.0 SUGGESTIONS FOR FUTURE WORK 


IITRI's suggestions for future work in this area entail 
1) refinements in the analytical techniques for unique materials, 
and 2) the assessment of space processing candidacy through con- 
trolled experiments and analysis. 

7 . 1 Improvements in Analytical/Empirical Techniques 

Improvements in our analytical-empirical prediction tech- 
niques can be derived from a continued cooperative effort between 
IITRI and North Americal Rockwell Corporation. Rockwell's part 
would be the laser melting of various oxide materials. IITRI 's 
part would be theoretical and analytical work leading to a quan- 
titative description of the glass forming ability of the oxides 
studied. 

R. A. Happe has indicated that in Rockwell's continuing 
program, "Manufacturing Unique Glar in Space" (UAS8-28991) , film- 
ing of the laser melting of six materials (V20^, Ge02, Nb20^, CoO, 
AI 2 O 2 , and ^20^) will be produced in a cooperative effort by Rockwell 
and Sandia Corporation. The purpose of these experiments will be 
to develop a densitometric film analysis technique for determining 
the temperature of the sample during the liquid-solid phase trans- 
formation. 

In IITRI' s analysis of the glass forming ability of .ne 
mullite system cr uhe current program, such a temperature-time 
history during the filmed laser melting experiment was not avail- 
able. Consequently, the required kinetic data had to be derived 
by either a) extrapolation, as for the viscosity-temperature rela- 
tion below the fusion temperature, or b) by deductive reasoning, 
as in the instance where we estimated the temperature where Happe 
observed the first detectable phase transformation in the mullite 
films. In the latter instance, "jiggling" of the cooling mullite 
sphere was observed during much of the crystallization process. 

The approximate viscosity that would be consistent with such a 
motion was estimated, thus leading to an estimated temperature 
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at which the event occurred. This analysis, although approximate, 
led to critical cooling rates for glass format-* on that are com- 
patible with empirical evidence (refer to Section 6.1). 

However, the cooperative Rockwell- Sandia effort to obtain 
an accurate temperature- time history during the laser melting 
experiments offers a \jnique opportunity for removing some of the 
assumptions and approximations in IITRI's analytical methods. 

For Instance, instead of assuming or extrapolating a shape for 
the viscosity-temperature relation of a given material below the 
liquidus, this relationship could be derived from analysis of the 
films. This analysis would proceed as follow : 

1. Starting with the temperature-coded films of the crystal- 
lization process, the entire experimentally observed growth 
rate-temperature curve would be generated. By counting 
frames, for instance, the growth over a small increment 

of time W’ould be computed. The calculated growth rate 
would then be referenced to the appropriate mean tempera- 
ture. This process would then be repeated at successive 
small time increments during the phase change process. An 
experimentally observed growth rate-ten?)erature relation 
would be generated. 

2. Assximing an analytical model for the growth process, 
the viscosity-temperature relation would then be gener- 
ated. This will enable us to obtain "glassy" data such 

as viscosity from materials that have never formed glasses. 

It is believed that the assumption of a particular growth 
model to give viscosity data in this analysis will give more 
accurate critical cooling rate data than the converse tech- 
nique of obtaining the growth behavior from an assumed 
viscosity model. Choice of which growth model ■. use might 
derive from morphological studies of crystallized samples 
as discussed in Sections 3.3. 1.3 and 3.3.2 above. 

3. Since growth and nucleation can be considered as independ- 
ent processes, the remainder of the analysis dealing with 
heterogeneous nucleation, homogeneous nucleation and critical 
cooling rates can proceed as in Section 6.1. 
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The materials for this proposed future work are listed In 
Table VII. Viscosity data for germanla, Ge02. are readily avail- 
able since Ge02 Is a well established ^^lass former. No viscosity 
data exists for the other materials. Therefore, the germanla 
viscosity data extracted from the literature will provide a good 
check on the accuracy of the above- described analytical-empirical 
technique to derive viscosity-temperature data below the fusion 
temperature. In addition, the viscosity- temperature data derived 
by this technique can be compared to the llquldus temperature 
viscosity data obtained by using an approximate technique des- 
cribed by Sokolov . 

7.2 Assessment of Space Processing Candidacy Through Controlled 

Experiments and Analysis 

To predict how a technique such as space processing could be 
utilized to produce glasses such as mulllte that are relatively 
unobtainable on Earth, a series c£ controlled laser melting ex- 
periments should now be conducted. The emphasis would be in deter- 
mining the optimum experimental conditions for obtaining near- 
intrinsic glass formation through in-space processing. Based on 
our Initial work and results described in Section 6.1, it is be- 
lieved that the mulllte system warrants this further work. 

The experimental variables for che laser melting experiments 
that TITRI suggests include- a) the amount of foreign material 
suc^ dust particles in the environment, and b) the amount of 
initial superheating neces' ary to dissolve siurfacc or internal 
refractory impurities into the melt without significant sample 
vaporization. Starting with films of the laser melting experiments 
conduct - under these conditions, IITRI would then perform analyses 
similar to those conducted on the current program. The result of 
this work will be a quantitative indication of how close to intrinsic 
the melting experiments were, and thus a quantitative estimate of the 
potential of space processing in producing unique glasses unobtain- 
able on earth. IITRI 's work on the current program serves as the 
basis for this prediction method. The Rockwell air suspended laser 
melting facility provides a unique capability for .ccomplishing this 
prediction of the space processing candidacy of mullite, and defining 

the experimental conditions that must be satisfied for the mission to 
be successful. 








8.0 CONCLUSIONS 



i 
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Prediction of the glass forming tendencies of unique oxide 
systems will facilitate determination of candidate materials for 
NASA's Space Processing Program. The ultimate objective of 
space processing is to produce technically significant glasses 
by either 1) extending the earth- limited regions of glass forma- 
tion f' r cert.'. r compositions, or 2) by obtaining glass formation 
in other compositions that are not glass formers based on empirical 
earth observations. 

IITRI has shown that these goals can be accomplished by 
moLtfiling the glass forming behavior of selected materials systems 
through the development of kinetic relationships describing nuclea- 
tion and crystallization phenomena. The theoretical/analytical 
background for this has been obtained on IITRI 's current program. 

A fairly comprehensive treatment of nucleation and growth kinetics 
in pure substances has been presented. The derived transformation 
kinetics have been successfully applied to a well-characterized 
system (SiO^) , an excellent glass former (820^) .. and a poor glass 
former by conventional means (Al20^) . The kinetic and thermodynamic 
parameters of viscosity and entropy of fusion were shown to be the 
primary materials parameters controlling glass forming cenden y. 

For complex multicomponent systems where diffusienal effects 
predominate, the state-of-the-art is not nearly as far a'^'''?ncecl as 
for simple substances. The transformation kinetics of mater ,.als 
which crystallize with a large compositional change are most prob- 
ably governed by the long range diffusion of the : '' zest- moving 
species. The general lack of specific diffusion data, however, 
dictate that simplifying assumptions be made, such as the validity 
of the Stokes-Einstein equation relating bulk diffusion and vis- 
cosity. 

In addition to diffusion-controlled kinetics, heterogeneous 
nucleation effects must be considered in real, non-simple systems, 
since heterogeneous nucleation is the probable cause of the earth- 
limited glass forming ability in many materials. Where complete 
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kinetic property data were not available for this analysis, we 
have shown that the required information can be derived from 
empirical evidence of the phase change phenomena. This empirical 
phase transformation information was obtained through a coopera- 
tive effort with North American Rockwell Corporation. The Rockwell 
experiments entailed the air-suspended laser melting of various 
oxide systems. The nature of the liquid-solid phase transformation 
(and thus glass forming tendency) was obtained from the analysis 
of the filmed records of these experiments . Using this analytical/ 
empirical approach the mullite system was analyzed, giving results 
that 1) are consistent with experimentally observed data, and 2) 
that indicated the promise of mullite as a future space processing 
candidate . 


This program provides NASA with a means of assessing the 
candidacy of selected systems for space processing, and determin- 


ing the experimental-hardware requirements of a successful in-space 
processing experiment. 


Respectfully submitted, 
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Research Engineer 
Ceramics Research Section 

APPROVED: 


/I 



S. A. Bortz.^, 

Assistant Director 

Mechanics of Materials Research Division 


92 





1 


REFERENCES 


1 . 

2 . 

3. 

4. 

5. 


6 . 

7. 

8 . 
9. 

10 . 

11 . 

12 . 

13. 


Morey, G, W. , The Properties of Glass, Reinhold. New York, 
2nd Edition (153^7: 

Jackson, K. A., "Nucleation from the Melt," Nucleation 
Phenomena , Americai Chemical Society (1966) . 

Turnbull, D. , "Thermodynamics and Kinetics of Formation of 
the Glass State and Initial Devitrification," Physics of 
Non-Crystall ine Solids, Proc. Int'l. Conf . , Delft, J. A. 
Prins. , ed. (l964) . 


Hammel, J. J., "Nucleation in Glass Forming Materials," 
r in N 
(19697 


Chapter in Nucleation, A. C. Zettlemoyer, ed. , Marcel D«Vker, 
N. Y. 


Hillig, W. B. , "A Theoretical and Experimental Investigation 
of Nucleation Leading to Uniform Crystallization of Glass," 
Symposium on Nucleation and Crystallization in Glasses and 
Melts, American Ceramic Society (1962). 

Collins, F. C., Z. Elektrochem. ^ 404 (1955). 

Turnbull, D. and Cohen, M. H. , "Crystallization Kinetics 
and Glass Formation," Modeini Aspects of the Vitreous State , 
Vol. 1, J. D. Mackenzie, ed. , But terworths , London ( 1960 ) . 

Uhlmann, D. R. and Chalmers, B. , "The Energetics of Nuclea- 
tion," Nucleation Phenomena , Am. Chem. Soc. (1966). 

Hammel, J. J. , "Nucleation in Glass - A Review," Adv. 
Nucleation and Crystallization in Glasses, Glass Division 
Symposium, Am. Ceram. Soc., Sp. Publ. #5, L. L. Hench and 
S. W. Freiman, eds . (1971). 

Rawson, H. , Inorganic Glass-Forming Systems , Academic Press, 
New York (1967) . 

Hillig, W. B. , "Glass as a Medium for Controlling Physio- 
chemical Reactions," in Reactivity of Solids . J. W. Mitchell, 
et.al. , eds., Wiley (1969). 

Uhlmann, D. R. , "Crystal Growth in Glass-Forming Systems - A 
Review," Adv. Nucleation and Crystallization in Glasses, 

Am. Ceram. Soc., Sp. Publ. #5, L. L. Hench and S. W. Freiman, 
eds. (1971). 

Uhlmann, D. R. , "A Kinetic Treatment of Glass Formation," 

J. Non-Cryst. Solids 7^ 337-48 (1972). 


93 




- “'■1 »• ! 



REFERENCES (cont’d) 


14. Mackenzie, J. D., J. Am. Ceram. Soc. ^(12) (1961) from Parks 
and Spaght, Physics 6(2) (1935) . 

15. IITRI Measured Container Composition; 71Si0.j • 15Na«0- 8CaO- 

4.5M 0-1.5AloO.,. ^ ^ 

g 2 3 

16. Reibling, E. F., J. Chem. Phys. ^(7) (1963) . 

17. "Recommended Values of the Thermophysical Properties of Eight 
Alloys, Major Constituents and Their Oxides," unpublished 
ccmpilation of data from Thermophysical Properties Research 
Center, Purdue University (1966); Source of Data: Kozakevitch, 
P., Met. Soc. Conf. 7 97 (1961). 

18. Fontana, E. H. and Pliimmer , W. A. , "A Study of Viscosity- 

Temperature Relationships in the G?.0., and SiO., Systems," Phys. 
Chem. Glasses 7 (4) (1966) . ^ ^ 

19. Bacon. J. F., et. al., "Viscosity and Density of Molten Silica 
and High Silica Content Glasses," Phys. Chem. Glasses 1 (3'^ 
(1960). 


20. Bruckner, R. , "Propertie.s and Structure of Vitreous Silica," 
J. Non-Cryst. Solids 5(3) (1970) . 


21. JANAF Thermochemical Data Tables. 

22. Bruckner, V. R. , Glass Tech. Ber. 37(H9) 413 (1964). 



k 


L 


I 

i 

I' 


E 

I 

I 

I 

I 


23. Macedo, P. B. , and Napolitano, A., J. Chem. Phys. 49(4) 1887 

(1968) . ~ 

24. Napolitano, A., Macedo, P. B. and Hawkins, E. G. , "Viscosity 
and Density of Boron Trioxide," J. Am. Ceram. Soc. 48(12) 

613 (1965). ~ 

25. Turnbull, D. and Cohen, M. H. , "Concerning Reconstructive 

Transformation and Formation of Glass," J. Chem. Phys, 29(5) 
(1958). ~ 

26. Hillig, W. B. , private communication. 

27. Hammel, J. J. , private communication. 

28. Frank, F. C., Proc. Roy. Soc. (London) A201 586 (1950). 

29. Takamori, T. and Roy, R. , "Rapid Crystallization of SiO^-Al.,0.5 

Glasses," J. Am. Ceram. Soc. ^(12) 639-44 (1973). ^ ^ ^ 


94 


i 





J 


I 


1 


1 





I 

I 

1 

1 

1 

I 

I 

I 

I 



I 

I 

V 

I 


i 

I 

I 


I 




REFERENCES (cont'd) 

30. Happe, R. A., "Manufacturing Uniaue Glasses in Space", 
Rockwell International Monthly Progress Report #11 to 
NASA-MSFC under Contract NASS -28991 (7 March to 5 April 
1974). 

31. Touloxikian, Y. S., "Recommended Values of the Thermo- 
physical Properties of Eight Alloys, Major Constituents 
and Their Oxides", unpublished data compilation at the 
Thermophysical Property Research Center (TPRC, Purdue 
Univ.); Data Source: Kozakevitch, P. Met. Soc, Conf. 

7 97 (1961). 

32. Brewer, L. , Professor of Materials Science, Univ. of 
Calif, at Berkeley, private communication. 

33. "Selected Values of Critical Thermodynamic Properties", 

NBS Circular 500. 

34. Happe, R. A., "Manufacturing Unique Glasses in Space", 
Monthly Progress Report No. 6 for the period of 1 October 
to 31 October 1973. NASA Contract NAS8-28991 (1973). 

35. ibid. Monthly Progress Report No. 8 for the period of 
1 December to 31 December (1973). 

36. Sargeant, P. T., "New Experimental and Theoretical 
Approaches to the Glassy State", PH.D. Thesis, Solid 
State Science, Pennsylvania State Univ., (Dec., 1967). 

37. Happe, R. A., "Study of the Production of Unique New 
Glasses", Final Report on Contract NAS8-28014 (estimated 
from information presented on pages 51-56) (13 June 1972) . 

38. Rawson, H., Inorganic Glass-Forming Systems , Academic 
Press, New York (1967) . 

39. Happe, R. A., "Manufacturing Unique Glasses in Space", 
NASA-MSFC Contract NAS 8-28991, Sixteenth Monthly Progress 
Report (17 Sept. 1974). 

40. Happe, R. A., "Manufacturing Unique Glasses in Space", 
NASA-MSFC Contract NAS8-28991, Seventeenth Monthly Progress 
Report (15 Oct. 1974). 

41. Ryshkewitch, E., Oxide Ceramics - Physical Chemistry and 
Technology , Academic Press, New York (i960) . 

4it. Happe, R. A., Rockwell International, private communication. 

43. Handbook of Chemistry and Physics , Chemical Rubber Co., 

46th Edition, pages D35-D37 (1964) . 

44. Sokolov, 0. K. , "Calculation of Viscosity in Molten Salts 
(Oxides)", NASA TTF-416 (circa 1962). 


95 





1 


1 


I 


4 


I 

I 

I 

I 

I 

I 


BIBLIOGRAPHY 


45. Avrami, M. , "Kinetics of Phase Change - II: Transformation- 
Time Relations for Random Distribution of Nuclei." J. Chem. 
Phys. 8 212 (1940). 

46. Berezhonoi, A. I., Glass Ceramics and Photo-Sitalls , Chapter 
II, "The Formation of Nuclei and the Crystallization of 
Glass," Mercol, S. A., Pincus, A. G. (trans.) Plenxim, New 
York (1970). 


T 

«> 


I 

I 

I 

I 

i 

I 

I 

! 

I 

I 

I 


47. Binsbergen, F. L. , "Computer Simulation of Nucleation and 
Crystal Growth," J. Cryst. Growth. 1^ 249 (1972). 

48. Binsbergen, F. L. , "Study of Molecular Phenomena in the Nuclea- 
tion of Crystallization by Computer Simulation and Application 
to Polymer Crystallization Theory," J. Cryst. Growth 13/14. 44 
(1972) . 

49. Buckle, E. R. , "Studies on the Freezing of Pure Liquids, II. 

The Kinetics of Homogeneous Nucleation in Supercooled Liquids," 
Proc. Roy. Soc. (London) A261 189 (1961). 

50. Calvert, P. D. and Uhlmann, D. R. , "Surface Nucleation Growth 
Theory for the Large and Small Crystal Cases and the Sip ifi- 
cance of Transient Nucleation," J. Cryst. Growth ^ 291 (1972). 

51. Cooper, A. R. , "Crystal Growth in Network Liquids by Structure 
Rearrangement," in Adv. Nucleation and Crystallization in 
Glasses, L. L. Hench and S. W. Frienan, eds.. Am. Ceram. Soc. 
Sp. Publ. #5 (1971). 

52. Cooper, A. R. and Gupta, P. K. , "Analysis of Diffusion 
Controlled Crystal Growth in Multicomponent Systems," ibid, 
pg. 131. 

53. Dunning, W. J. , "General and Theoretical Introduction (to 
Nucleation," Chapter in Nucleation , A. C. Zettlemoyer, ed. , 
Marcel Dekker, N. Y. (1969) . 

54. Eckstein, B. , "Vitreous States," Mat. Res. Bull. 3 199-208 
(1968) . 

55. Ewing, R. H. , "The Free Energy of the Crystal-Melt Interface 

from the Radial Distribution Function," J. Cryst. Growth 11 
221 (1971). ~ 

56. Fisher, J. C., et. al. , "Nucleation," J. Appl. Phys. 775 
(August 1948) . 


I 

I 


96 



BIBLIOGRAPHY (Cont ' d . ) 


57. Frank, F. C. , "An Outline of Nucleation Theory," J. Cryst. 
Growth 13/14 154 (1972) . 

58. Ghez, R. , "An Exact Calculation of Crystal Growth Rates Under 
Conditions of Constant Cooling Rate," J. Cryst. Growth ^ 

153 (1973). 

59. Ghez, R. and Lew, J. S., "Interface Kinetics and Crystal 
Growth Under Cor.ditons of Constant Cooling Rate," J Cryst. 
Growth ^ 273 (1973) . 

60. Gibbs, J. H. , "Nature of the Glass Transition and the Vitreous 
State," in Modern Aspects of the Vitreous State, I , J. D. 
Mackenzie, ed. 

61. Gilmer, G. H. and Bennema, P. , "Computer Simulation of 
Crystal Surface Structure and Growth Kinetics," J. Cryst. 

Growth 13/14 148 (1972) . 

62. Gutzow, I. and Kashchiev, D. , "Kinetics of Overall Crystal- 
lization of Undercooled Melts in Terms of the Non-Steady 
State Theory of Nucleation," in Adv. Nuclation and Crystal - 
lization in Glasses , L. L. Hench and S. W. Freiman, eds.. Am. 
Ceram. Soc. Sp. Publ. #5 (1971). 

63. Hammel, J. J. , "Direct Measurements of Homogeneous Nucleation 
Rates in a Glass-Forming Systems," J. Chem. Phys. ^ (6) 

2234 (1966). 

64. Heaton, H. M. and Moore, H. , "A Study of Glasses Consisting 
Mainly of the Oxides of Elements of High Atomic Weight, 

Part III. The Factors Which Determine the Possibility of 
Glass Formation," J. Soc. Glass Technol. XLI (198) (Feb. 1957). 

65. Hillig, W. B. , "Kinetics of Solidification from Nonmetallic 
Liquids," in Kinetics of High Temperature Processes , W. D. 
Kingery , ed . 

66. Hillig, W. B. and Turnbull, D. , "Theory of Crystal Growth in 
Undercooled Pure Liquids," J. Chem. Phys. (4) 914 (1956). 

67. Hopper, R. W. , et. al. , "Crystallization Statistics, Thermal 
History and Glass Formation," J. Non-Cryst. Solids ^5-62 
(1974). 

68. Hopper, R. W. and Uhlmann, D. R. , "Solute Redistribution 
During Crystallization at Constant Velocity and Constant 
Temperature," J. Cryst. Growth ^ 203 (1974). 


97 




BIBLIOGRAPHY (Cont'd.) 



69. Hopper, R. W. and Uhlmann, D. R. , Temperature Distributions 
During Crystallization at Constant Velocity," J. Cryst. 

Growth 19 111 (1973). 

70. Jackson, K. A., "On Surface Free Energy Calculations," J. 
Cryst. Gro/^th ^ 119 (1971). 

71. Jones, D. R. H. , "A Mechanical Model of Solid-Liquid Inter- 
faces," J. Cryst. Growth ^ 167 (1974). 

72. Knight, C. A. , "Breeding of Crystal Nuclei by Classical 
Nucleation: Theory and Some Observations and Experiments," 

J. Cryst. Growth U 201 (1971). 

73. Konak, A. R. , "Difficulties Associated with Theories of Growth 
from Solution," J. Cryst. Growth 247 (1973). 

74. Lothe, J. and Pomd, G. M. , "Statistical Mechanics of Nuclea- 
tion," Chapter in Nucleation, A. C. Zettlemoyer, ed. , Marcel 
Dekker, N. Y. (19697: 

75. MacDowell, J. F. , "Nucleation in Glasses," in Nucleation 
Phenomena . Am. Chem. Soc. (1966). 

76. Magill, J. H. and Li, H. M, , "A Corresponding States Equation 
for Crystallization Kinetics," J. Cryst. Growth 361 (1973). 

77. Markov, I. and Kashchiev, D. , "Nucleation on Active Centres 1. 
General Theory," J. Cryst. Growth 1^ 170 (1972). 

78. McMillan, P. W. , Glass-Ceramics , Chapter 2, "Crystallization 
and Devitrification," Academic Press, N. Y. (1964). 

79. Robertson, D. and Pound, G. M. , "Numerical Simulation of 

Heterogeneous Nucleation and Growth, J, Cryst. Growth 19 269 
(1973). ^ 

80. Sargeant, P. and Roy, R. , "New Approach to the Prediction of 
Glass Formation," Mat. Res. Bull. 3 265-80 (1968). 

81. Sears, G. W. , "Recent Developments in Nucleation Theory," in 

, C. Klingsberg, ed. , Gordon 


82. Stevels, J. M. , "Permitted and Forbidden Structural Units in 
Vitreous Systems," Mat. Res. Bull. 3 599-610 (1968). 

83. Sun, K. , "Fundamental Condition of Glass Formation," J. Am. 
Ceram. Soc. 30 (9) 277 (1947). 


Physics and Chemistry 
and Breach (1963) . 


of Ceramics 


98 





BIBLIOGRAPHY (Cont'd.) 


84. Tashiro, M. , "Nucleation and Crystal Growth in Glasses," 
Eighth International Congress on Glass (1968) . 

Turnbull , D. , "Phase Changes," in Solid State Physics - 

Advance.^ in Research and Applications, Vol . 3 . Academic Press, 


86. Turnbull, D., "Under What Conditions Can a Glass be Formed?" 
Contemp. Phys. 10 (5) 473-88 (1969). 

87. Turnbull, D. and Cohen, M. H. , "Free-Volume Model of the 
Amorphous Phase: Glass Transition." J. Chem. Phys. 34 (1) 

120 (1961). “ 

88. Walton, A. G. , "Nucleation in Liquids and Solutions," Chapter 
in Nucleation, A. C. Zettlemoyer, ed. . Marcel Dekker, N. Y. 

(i9^?y: 

89. Walton, A. G. , "Bulk and Surface Transport in Liquids - An 
Introudctory Treatment," Contemp. Phys. (5) 489 (1969). 

90. White, W. B. and Berkes , J. S., "Calculation of Volume Free 

Energy of Nucleation of a Crystalline Phase from a Multi- 
component Regular Solution," J. Am. Ceram. Soc. 52 (4) 231 
(1969). ~ 


•I 






99 



